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Abstract. Experimental progresses in the miniaturisation of electronic devices have 
made routinely available in the laboratory small electronic systems, on the micron 
or sub-micron scale, which at low temperature are sufficiently well isolated from 
their environment to be considered as fully coherent. Some of their most important 
properties are dominated by the interaction between electrons. Understanding their 
behaviour therefore requires a description of the interplay between interference effects 
and interactions. 

The goal of this review is to address this relatively broad issue, and more specifically 
to address it from the perspective of the quantum chaos community. I will therefore 
present some of the concepts developed in the field of quantum chaos which have some 
application to study many-body effects in mesoscopic and nanoscopic systems. Their 
implementation is illustrated on a few examples of experimental relevance such as 
persistent currents, mesoscopic fluctuations of Kondo properties or Coulomb blockade. 
I will furthermore try to bring out, from the various physical illustrations, some of the 
specific advantages on more general grounds of the quantum chaos based approach. 
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1. Introduction 

The title of this review may seem self contradictory in two respects. To begin with, 
it associates chaos, a purely classical notion, with quantum physics. Furthermore it 
implies that this association, which as we will see refers traditionally to the study of 
low-dimensional non-interacting quantum systems, will be considered in the context of 
many-body physics. 

The first of these contradictions is however mainly a question of semantics. Indeed, 
if in an early period of development of the field of quantum chaos, some of the issues 
addressed had to do with the possible existence of true chaotic dynamics for quantum 
systems, this was relatively quickly answered, essentially in the negative. Quantum 
chaos now mainly refers to the study of the consequences, for a quantum system, of the 
more or less chaotic nature of the dynamics of its classical analogue. It has followed two 
main avenues. The first one is based on semiclassical techniques - specifically the use 
of semiclassical Green's function in the spirit of Gutzwiller's trace formulae (Gutzwiller 
1971, Balian & Bloch 1972, Berry & Tabor 1977, Gutzwiller 1991, Gutzwiller 1990) - 
which provides a link between a quantum system and its 7i — > limit. The second is 
associated with the Bohigas Giannoni Schmit conjecture (Bohigas et al. 1984a, Bohigas 
et al. 19846, Bohigas 1991), or related approaches (Peres 1984), which states that the 
spectral fluctuations of classically chaotic systems can be described using the proper 
ensembles of random matrices. 

Some of the beauty of quantum chaos is that it has developed a set of tools 
which have found applications in a large variety of different physical contexts, ranging 
from molecular and atomic physics (Wintgen & Friedrich 1986, Delande & Gay 1986, 
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Wintgen 1987, Delande et al. 1991), to acoustics (Derode et al. 1995, Fink 1997, Fink 
et al. 2000), nuclear physics (Bohigas & Leboeuf 2002, Olofsson et al. 2006), cold atoms 
(Mouchet et al. 2001, Hensinger et al. 2001, Steck et al. 2001, Steck et al. 2002), optical 
(Nockel & Stone 1997, Gmachl et al. 1998) or microwave (Stockmann & Stein 1990, 
Sridhar 1991, Kudrolli et al. 1995, Alt et al. 1995, Pradhan & Sridhar 2000) resonators, 
and of course mesoscopic physics (Richter et al. 19966, Richter 2000, Alhassid 2000). 
With few exceptions (see nevertheless (Bohigas & Leboeuf 2002, Olofsson et al. 2006)) 
most of these physical systems share the property of being correctly described by non- 
interacting, low- dimensional, models. 

This is true in spite of the fact that random matrix ensembles were introduced in 
the early fifties by Wigner (see the series of articles reprinted in Porter 1965) to explain 
the statistics of slow neutron resonances, and were therefore applied in the context of 
strongly interacting systems. In that case however, it was less the notion of chaos than 
the one of complexity (large number of degrees of freedom, strong and complicated 
interactions) which was proposed by Wigner to justify this approach. 

At any rate, the scope of this review will be concerned with a very different 
type of interacting system, namely the Landau-Fermi liquid, for which the system 
can be explain as a set of quasi-particles interacting weakly through a (renormalized) 
interaction amenable to perturbative treatment. More specifically, what we have in 
mind are various realizations of fully coherent, confined electron gasses, with a density 
high enough that a Landau-Fermi-liquid type description applies. These are typically 
semiconductor quantum dots or small metallic nano-particles, within which a few tens 
to a few thousands of electrons interact through a screened Coulomb interaction. 

Although this screened interactions between electrons is weak, and is therefore well 
described by a standard perturbative approach (first order perturbation in the simplest 
cases, or eventually with some re-summation of higher order terms in other situations), 
some important physical processes are actually largely dominated by them. Moreover, 
the systems considered are only weakly affected by their environment and can therefore 
be assumed fully coherent. Because of the confinement, translational invariance is then 
broken, and some new and interesting physics is brought in by the fact that, in the 
non-interacting limit from which the perturbation scheme is developed, eigenstates are 
not just plane waves. The mesoscopic fluctuations associated with confinement and 
interference need to be taken into account for the eigenstates and one-particle energies, 
either at a statistical level or in a detailed way associated with a given geometry. 

How to describe these mesoscopic fluctuations, and implement their consequences 
for many-body effects can be done in a variety of ways. For diffusive systems, techniques 
based on diagrammatic perturbation expansion in the disorder potential can be used 
(Altshuler & Aronov 1985, Aleiner et al. 2002, Akkermans & Montambaux 2007), as well 
as approaches based on the supersymmetric cr-model (Efetov 1999, Mirlin 2000), which 
is also appropriate for the description of ballistic chaotic systems (Blanter et al. 2001a) 
(see also (Andreev et al. 1996, Leyvraz & Seligman 1997, Agam et al. 1997) in this 
context). In this review, I shall however limit myself to the methods coming from 
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the quantum chaos community. One reason for this limitation of scope is that there 
already are very good and complete reviews which give a excellent account of the other 
approaches. Another is that the quantum chaos perspective is in many useful cases 
more intuitive, and somewhat simpler to apply from a technical point of view. As a 
consequence, this will make it possible to present in an essentially self-contained way 
the technical tools required to understand a large class of many-body effects relevant for 
these quantum dots or nano-particles, as long as they are in the Landau- Fermi-liquid 
regime. The goal is that it should be possible to follow almost all of the review with a 
graduate level in quantum mechanics and, in some cases, basic notions of many-body 
theory such as what can be found in classic textbooks such as Fetter and Walecka (Fetter 
& Walecka 1971). This, I hope, will make it particularly convenient for experimentalists 
or theoreticians who wish to enter into this field. 

Another advantage of the quantum-chaos based approach is that it is by nature 
more flexible, and is therefore not limited to chaotic or diffusive dynamics. How much 
physics is missed by other points of view because of this limitation can be debated, and 
I will return to this discussion at the end of this review. However, if for metallic nano- 
particles the choice of a description in term of disordered (diffusive) system is dictated 
by the actual physics of these materials, it is clear that for semiconductor quantum dots, 
one reason for why so much focus has been put on chaotic dynamics is that it is the only 
one that could be addressed by the more traditional techniques of solid state physics. 
Having a tool which makes it possible to consider other kinds of dynamics at least gives 
the possibility of asking the question of whether anything new, or interesting, can be 
found in these other regimes. 

The structure of this review is therefore the following. A first section will be devoted 
to the description of the basic tools necessary to study the physical problems we want 
to consider. As we want to address different energy scales, or from an experimental 
point of view different temperature ranges, it will be necessary to introduce a few 
complementary point of views. Semiclassical techniques, and in particular the use of 
semiclassical Green's function, will be well adapted to temperature ranges significantly 
larger than the mean level spacing A. They will be the subject of section 2.1. The 
low (T < A) temperature regime however requires a modeling of individual energies 
and wave-functions, and are therefore better described, in the hard chaos regime, by 
statistical approaches such as random matrix theory and the random plane wave model. 
These latter will be introduced in section 2.2. Finally section 2.3 will provide a discussion 
of the screening of the Coulomb interaction. 

I will then turn to the description of a few examples of physical systems where the 
physics is dominated by the interplay of interaction effects and mesoscopic fluctuations. 
The choice of these examples is of course quite arbitrary, and the criterion I have used 
to select them is essentially my familiarity with the issue. Therefore, there will be a 
strong bias toward questions I have actually worked on, which should not be interpreted 
as a statement about their relative importance. I will start with a discussion of the 
orbital magnetic response, with some general considerations in section 3 followed by 
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a few specific examples of diffusive and ballistic systems in section 4. One important 
difficulty to be addressed here is the renormalization of the interaction due to higher 
order terms in the Cooper channel, and this issue will be discussed in detailed in both 
the diffusive and ballistic regimes. In section 5, 1 will wander briefly away from Fermi 
liquid systems and address the mesoscopic fluctuations associated with the physics of 
a Kondo impurity (Kondo 1964, Hewson 1993) placed in a finite-size system. The last 
physical example will be, in Section 6, the role of interactions in the fluctuation of 
peak spacing in Coulomb blockade (Beenakker 1991, Weinmann et al. 1996, Grabert 
& H. Devoret 1992, Sohn et al. 1997) experiments. After a general introduction of 
the universal Hamiltonian picture, I will cover the various physical effects which needs 
to be further considered if one expect to understand experimental peak spacing and 
ground-state spin distributions for these systems. 

Finally, the concluding section will contain some general discussion, and in 
particular will come back to the issue of non-chaotic dynamics. 

2. Basic tools 

2. 1 . Semiclassical formalism 

Consider a system of indistinguishable Fermions governed by the one-particle 
Hamiltonian 

and interacting weakly through the two body potential V(r — r'). A systematic 
perturbative expansion can be constructed to arbitrary order (if necessary) in terms 
of the unperturbed Green's function f 

G(ry ;e )*(r| -L^ M = £*^MM, (2 . 2) 

where in the last expression e K and (p K are respectively the energies and eigenstates of 
-f/ip. In a clean bulk system U(r) = so that the eigenstates are just plane waves, and 
the expression of the Green's function becomes trivial in the momentum representation. 
For confined (coherent) systems however, translational invariance is lost and there is in 
general no simple expression for the exact eigenstates and eigenfunctions. It therefore 
becomes necessary to find some approximation scheme for the unperturbed Green's 
function itself before considering a perturbation expansion in the interaction. 

2.1.1. Semiclassical Green's function In many regimes of interest a semiclassical 
approach can be used to fulfil this role. This includes the case where the confining 
potential U(r) is a smooth function on the scale of the Fermi wavelength Af, but also 
for instance when it contains a weak, eventually short range, disorder, as long as Af is 

| All the Green's functions used in this review will be unperturbed Green's functions. I shall therefore 
not use any subscript to distinguish them from the interacting ones. 
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much smaller than the elastic mean free path £. Under these conditions, the retarded 
Green's function G R (e) == lim^o G(e + irj) can be written as a sum over all classical 
trajectories j joining r' to r at energy e (Gutzwiller 1990, Gutzwiller 1991) 

G R (r, r'; e) ~ £ Gf (r, r'; e) 



Gf(r,r' ;e ) 



def 



2tt 



(2i7T^)( d+1 )/ 2 

with d the number of degrees of freedom, 



p ■ dr 



the classical action along the trajectory j, and 



1 



■det' 



d 2 S 
drdr' 



1/2 



(2.3) 



(2.4) 



(2.5) 



a determinant describing the stability of trajectories near j. In (2.5) r = (r±, ■ ■ ■ ,r^), 
and the prime on the determinant indicates that the first component (i.e. first row 
d 2 S/dr 1 dr' i and first column d 2 S/dr i dr' l) % — 1, • • • , d) is omitted. Finally, the Maslov 
index Q essentially counts the number of caustics (i.e. places where the determinant Dj 
is zero) on the trajectory j between r' and r §. For two dimensional systems (d — 2), 
the determinant takes the particularly simple form 

V2 



(2.6) 



ryrjl (dr±/dp' ± ) r± 

where ry and rj_ are the r-components respectively parallel and orthogonal to the 
trajectory. 



2.1.2. Simple properties of the classical action Many important characteristic features 
of the semiclassical Green's function, and therefore of the fermion g be directly 

deduced from basic properties of the classical action (Arnold 1989, Goldstein 1964). In 
particular: 

i) the variation with respect to the energy 



de 



t 



3 i 



(2.7) 



where tj is the time elapsed to go from r' to r along trajectory j at energy e; 
ii) the variation with respect to the position 



dr' ^ dr 
and finally 

iii) the effect of a perturbation. 

§ For the kind of kinetic plus potential Hamiltonian we consider here, the Maslov index increase by 
one at each crossing of a caustic. Note however that for more general Hamiltonian it may however also 
decrease. 



(2.8) 
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Indeed, let us assume that the one particle Hamiltonian can be written as the sum of a 
main term H and a small perturbation Hi 

Hi p = H + Hi, (2.9) 

and let us denote by S® the action calculated for the trajectory j, i.e. (r°(£),p°(i))j, 
t E [0,t°], joining r' to r under the Hamiltonian H . We then have 

^ = Sj - S° ~ - f ° dtHiir^t)^)) . (2.10) 

Note there is no need for ffi to be small on the quantum scale, and therefore (2.10) 
remains applicable much beyond the limit of quantum perturbation theory. || 

To illustrate how the above properties can be used in our context, let us consider 
for instance the (unperturbed) local density of states, 

(r; e) d = f £ |y? K (r)| 2 5(e - e K ) = --ImG>, r; e) . (2.11) 



Using (2.3), v\ oc can be expressed as a sum over all the closed orbits starting and ending 
at the point r. In this process, the "direct" orbit jo, whose length goes to zero as 
r — > r', needs however to be treated separately, as the determinant D jo diverges. On 
the other hand, the contribution of this orbit to the Green's function can be identified 
to the free Green's function for a constant potential. It therefore just gives rise to the 
smooth (bulk-like) contribution 



/ up . , , m e k(r) d 2 d7r d/2 

■ w f r :e)= / , *M e-H(p,r))= , \ ' n -— -. (2.12) 

1 1 (2-)'7> 2 r(d/2 + l) 



/(2^F 5(e " 



(i/q 2 '' = m e /27rh 2 , Uq 3 ^ = m e k/2ir 2 h 2 ), where d is the dimensionality and the last equality 
holds for the usual kinetic plus potential Hamiltonian, with k(r) = i/2m e (e — U(r))/H. 
I shall in the following use the notation p(e) = f dru\ oc (r; e) for the total density of 
states, and 

Po{e) = I§W 5{e - H ^ r)) (2 ' 13) 
for its Weyl (smooth) part. 

The local density of states can therefore be separated into a smooth and an 
oscillating part 

v\ oc (r; e) = i/ (r; e) + v osc (r; e) (2.14) 
where z/ osc (r; e) is expressed semiclassically as a sum over all finite length closed orbits 
^osc(r-;e)= ^ Vj(r;e) + c.c. (2.15) 

C) = (2mh)(*+W D M 6XP ^ S ^l k ~ <^/ 2 ) • ( 2 - 16 ) 

|| It should be stressed that the exponential sensitivity to perturbations of chaotic trajectories does 
not prevent finding a perturbed trajectory following closely the unperturbed one and joining the same 
cndpoints in configuration space (but with a slightly different momenta). See for instance Cerruti & 
Tomsovic (2002) for a recent discussion of this (old) question. 
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For energies e close, on the classical scale, to some reference energy e, one can 
therefore use (2.7) to write 

Uj(r;e + 8e) = Uj(r;e) exp(i8etj/h) . (2-17) 

Thus, the local density of states appears as the bulk contribution plus some oscillating 
terms which, with (2.7), have a period in energy 2nh/tj determined by the travel time 
of the closed orbits. 

In the same way, Friedel oscillations near the boundary of the system or near an 
impurity can be understood as a direct consequence of (2.8), applied to the trajectory 
bouncing on the obstacle and coming back directly to its starting point. Quite generally 
one can write for the contribution of the orbit j to the local density of states 

Vj(r + 5r; e) = Vj(r; e) exp(i(p^. - Pj)5r/h) , (2.18) 

so that locally v\ oc {r;e) appears as a sum of plane waves the wave vectors of which 
are determined by the difference between the final and initial momentum of the 
corresponding returning orbits. Periodic orbits, which are such that p'j = pj, have 
no variation locally and therefore will give rise to the dominant contribution to the 
total density of states p(e) (Gutzwiller 1970, Gutzwiller 1971). Friedel oscillations on 
the other hand correspond to trajectories which, after bouncing off the boundary of the 
system or some impurity, travel back directly to the initial point, so that p'j = —Pj- In 
this case the corresponding plane wave contribution has a wave vector with modulus 
twice the wave vector k{r) = y/(e — U(r))/2m e /H. In the particular case of a two 
dimensional system with a straight hard wall (with Dirichlet boundary condition) at 
x — 0, direct application of (2.6) and (2.15) gives (for kx > 1) 

■Wr=(,, y );e) = -v^f 3in(2 ^ /4) , (2.19) 

from which the Friedel oscillations in the density of particles n(r) is derived by 
integration over the energy. More general cases (e.g. curved boundary) are easily 
obtained by calculating the corresponding value of dp±/dr' ± (and eventually Maslov 
indices). 

Finally as an illustration of the third property (2.10), let us compute the variation 
of the density of states under the modification of an external magnetic field B = V x A. 
When the magnetic field is changed from B — > B + 5B (with the corresponding change 
of the vector potential A — > A + 5 A) , one has 

H = ^-(p-eA) 2 + U(r), (2.20) 

and in first order in 5 A 

H x = —(p-eA)-5A = ev-5A, (2.21) 

m e 

with v = r = (dH/dp) the velocity. The variation of the action along a closed trajectory 
j : r — > r is therefore given by 

\SS 6 = \<[> dt(5A-v) = l<( dl-5A = ^^, (2.22) 
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where, 5<p j is the flux of 5B across the trajectory j and O = 2irh/e is the quantum 
flux. The variation of the contribution of the orbit j to the local density of states can 
therefore be written as 

Vjir; e;B + 5B) = Uj(r; e; B) exp(i<%/0 o ) • (2.23) 

2.1.3. Sum rule for the determinant Dj The computation of the contribution of some 
orbit j to the semiclassical Green's function for a given geometry implies in practice the 
determination of the action of the orbit, which is usually not too difficult, but also of the 
stability determinant Dj and the Maslov index Q which for generic systems may involve 
some technicalities (see e.g. Bogomolny (1988) for an illustration on the examples of 
the stadium and elliptic billiards, and Creagh et al. (1990) for a detailed discussion 
about the Maslov indices). It turns out that in practice a large number of results can be 
obtained without an explicit calculation of these quantities, but can be derived from a 
sum rule (M- formula) for the determinants Dj, analog in spirit to the Hannay - Ozorio 
de Almeida sum rule (Hannay & Ozorio de Almeida 1984, Ozorio de Almeida 1988), and 
which in a similar way expresses the conservation of Liouville measure by the classical 
flow. The M-formula can be expressed as (Argaman 1996) 

E T^W-V ="i? ) (r>)P; i (ry,t) , (2.24) 

where v^\r') is the bulk density of states per unit area (and spin) (see (2.12)) for the 
local value of k F and P^(r, r', t) is the classical (density of) probability that a trajectory 
launched in r' is in the neighbourhood of r at time t. 

The sum rule (2.24) is particularly useful for diffusive systems, for which P^(r, r', t) 
is solution of a diffusion equation (with diffusion coefficient D) 

(d t + DA r )P^(r, r', t) = 5{r - r')S(t) (2.25) 

with boundary conditions daP c \ = at the boundary of the system (if any). We shall see 
in the next subsection that it can be also applied usefully for ballistic chaotic systems. 

2.1.4- Thouless Energy When considering a confined system of (for now) non- 
interacting particles, one might first, before any actual calculation, try to understand 
what are the energy scales affected by the confinement. On the low-energy side, this is 
clearly bounded by the mean level spacing, the finiteness of which is the most obvious 
consequence of the fact that the system is bounded. On the high energy end, a direct 
implication of (2.17) is that if the Green's function is smoothed on an energy scale e, 
only the contributions of trajectories with a time of travel t < K/e survive the averaging 
process. Therefore the minimal time t m j n such that a classical particle feels the presence 
of the boundary determines the maximal energy scale E^h such that the quantum system 
is affected by this latter. This energy scale, £rh, is referred to as the Thouless energy. 

For ballistic systems, t m j n is essentially the time of flight t^ across the system, which 
is also the time scale of the shortest closed orbit for a generic point r inside the system. 
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As a consequence E Th is also in this case the energy scale beyond which no fluctuations 
exist. 

For diffusive systems t m i n is the time necessary to diffuse to the boundary of the 
system, time at which the solution of the diffusion equation (2.25) starts to differ from a 
free space diffusion. In that case however, this scale is different (and usually significantly 
larger) than the time associated with the shortest closed orbit, which is rather of the 
order of the momentum randomisation time t tl . 

As I shall illustrate below, a significant number of results can be derived for singly 
connected chaotic or diffusive quantum dots using the sum rule (2.24) in conjunction 
with a simple approximation for the probability P^(r,r',t). Indeed, for time shorter 
than t min the presence of the boundaries can be ignored and the free flight (free 
diffusion) expression can be used for the ballistic (diffusive) case. On the other hand 
the classical probability P^(r,r',t) can usually be taken as being independent of the 
initial condition for time larger than t min and simply proportional to the phase space 
volume J dp 5 (e — H(p, r)). For a billiard system of volume A this gives for instance 
for the return probability 

P^(r,r,t) = t<t min (2.26) 

P^r,r,t) = l/A t>t min (2.27) 

for chaotic system, and in the diffusive case 

P c e ,(r, r, t) = (AnDt)" 1 / 2 t < t min (2.28) 
P^(r,r,t) = l/A t>t min . (2.29) 

2.2. Random matrix and random-plane-wave models in the hard-chaos regime 

The semiclassical approach introduced in the previous subsection is the natural tool to 
describe energy scales significantly larger than the mean level spacing A. It is however 
not convenient, would it be only because it is usually not convergent, when the properties 
of a single wave-function are considered, and more generally when quantum properties 
are investigated on the scale not larger than a few mean level spacings (see however 
Tomsovic et al. (2008) in this respect). 

In the hard-chaos regime (and quite often in the diffusive regime), it is however 
possible to use an alternative route, based on the statistical description of the eigenstates 
and eigenf unctions fluctuations. 

2.2.1. Random matrices The most basic model is to assume that the fluctuations of 
physical quantities for the quantum system under consideration are well described by 
ensembles of random matrices, such as the Gaussian orthogonal, unitary or symplectic 
ensembles. These ensemble have been introduced by Wigner in the context of nuclear 
physics to account for the complexity (i.e.: large number or degrees of freedom, 
large and complicated interactions) characteristic of the nuclei. Studies of billiard 
systems have however led Bohigas, Giannoni and Schmit (Bohigas et al. 1984a, Bohigas 
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et al. 1984 b) to conjecture that even "simple" (i.e.: low dimensional, with innocent 
looking Hamiltonians) systems would display the spectral fluctuations of these Gaussian 
ensembles as long as they have a chaotic dynamics. This conjecture, although still 
not formally proven, is supported by recent semiclassical calculations showing that the 
form factor (the Fourier transform of the two point correlation function) predicted 
by the random matrix ensembles can be recovered in all order of a perturbation 
expansion within a periodic orbit theory (Berry 1985, Richter & Sieber 2002, Muller 
et al. 2004, Stefan Heusler & Haake 2004, Muller et al. 2005). I has furthermore 
been verified numerically on a large variety of systems, giving a kind of experimental 
demonstration that classical chaos, rather than complexity, is the origin of these 
characteristic fluctuations properties. 

Wigner Gaussian ensembles are constructed by first considering the set of N x N 
(in the limit N — > oo) Hamiltonian matrices corresponding to the symmetries with 
respect to time reversal of the physical systems. Those are Hermitian matrices H = 
H^+iH® for time-reversal non-invariant systems (Gaussian unitary ensemble (GUE)), 
real symmetric matrices H = for spinless time-reversal invariant systems (Gaussian 
orthogonal ensemble, (GOE)) , and quaternion real matrices H = <g> 1 + <g> a x + 

spin 1/2 time- reversal invariant but non-rotationally invariant 
systems (Gaussian symplectic ensemble, (GSE)). (H^ is a real symmetric matrix, and 
the H^ a \a > are real antisymmetric matrices.) The associated probability is then 
constructed i) assuming that each matrix element hff (i < j) is independent; ii) in such 
a way that the probability is invariant under the group transformation corresponding 
to a change of basis (unitary transformations for GUE, orthogonal transformations for 
GOE, and symplectic transformations for GSE). This leads to (Mehta 1991) 

P PnMT (H)dH = K Nf3nMT exp [-Tr(H 2 )/4v 2 ] dH (2.30) 

where -fGv/3 RMT * s a normalization constant, /3rmt indexes the symmetry class (/3rmt — 
1,2,4 corresponding respectively to GOE, GUE and GSE), v is an energy scale 
determined in practice by the physical mean level spacing, and dH is the natural measure 

i<j a 

From the probability distribution (2.30) various spectral correlation functions can be 
derived (Mehta 1991). For instance, the distribution of the (scaled) nearest neighbor 
s = (e n+ i — e n )/A can be shown to be well approximated by the Wigner-surmise 
distribution 

Pnns(s) = a &MT /- T exp(- C/3RMTS 2 ) , (2.32) 

with the numbers (ag RMX , Cp RMT ) fixed by the constraints on the normalization and on 
the mean ^f. 

The parameters (a/3 RMT , %„ T ) are equal to (7r/2,7r/4) for /3 RM t = 1, (32/tt 2 , 4/tt) for /3 RM t = 2, 
and (2 18 /3 6 7r 3 ,64/97r) for /3 RMT = 4 (see e.g. Appendix A of Bohigas (1991)). 
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The main content of (2.30) is its universal character. Indeed beyond the scale v and 
the symmetries of the system, the resulting distributions are completely independent 
of the particular feature of the physical problem under consideration, as long as 
the corresponding classical dynamics is chaotic. This makes it possible to obtain 
quantitative predictions for various physical configuration without a precise knowledge 
of the system specific details. This is presumably one of the reasons why so much focus 
has been put, both theoretically and experimentally, on the hard chaotic regime. 

2.2.2. Random-plane-wave model As long as spectral statistics are concerned, and 
that energy scales much shorter than the Thouless energy are considered, the random 
matrix models have been shown to be extremely reliable. The situation is however more 
ambiguous when one considers wave-functions statistics. On one hand, some properties, 
like the Porter-Thomas character (Brody et al. 1981) 

P(u= \<p n (r)\ 2 /(\<p\ 2 )) = exp(-u) GOE (2.33) 

= ^^exp(-w/2) GUE (2.34) 

V 27TM 

of the fluctuations around the mean value (|y?| 2 ) of the eigenstates probability at 
a given position, are well observed in numerical calculations, and can be derived 
straightforwardly from a random matrix description. On the other hand some 
other wave-function statistics clearly cannot be addressed within a simple random 
matrix model. Consider for instance a billiard-like quantum dot, which is therefore 
characterized by a constant (in space) Fermi wave vector k F . Wave-functions correlations 
are then characterized by a length scale Af = 2n/k F which is clearly absent in the 
random matrix description. 

To introduce this scale in the wave-function statistics let us consider, for an 
arbitrary wave function (p(r), its Wigner transform defined as 

[(p] w (r,p) d = J dxexp( y —ipx/h)(p*( y r + x/2)(p(r — x/2). (2.35) 

The normalization of the wave function implies (27rH)~ d J drdp[<p n ]-w = 1. If tp is an 
eigenstate with energy e, it can be shown (Berry 1977, Voros 1979) that for chaotic 
systems [<p\w( r ,P) converges in probability in the semiclassical limit toward the micro- 
canonical distribution 

AW(r,p) = p 1 5(e n - H(r,p)) (2.36) 

(p is the Weyl density of states (2.13).) the definition (2.35) can be inverted into 

V*( r W> = (2^Y* / tpexPbP ■ (r - 0/fi)Mw(r,p) ■ (2.37) 

Replacing on average [</?]w by /i mc one immediately obtains an explicit expression for 
the two-point correlation function (ip*(r)<p(r')). For instance for the usual kinetic plus 
potential Hamiltonian H = p 2 /2m + U(r), and for distance short enough that the 
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variation of the potential (and therefore of k(r)) can be neglected 

i *( \ ( >\\ u ^ r ) d0 ci\ 'i (aw 
(ip [r)ip(r)) — / — exp(ik\r — r \ cos{0)) 

ft Jo ^ 



<a 4 Xv> = ^"W • (2-40) 



= *^-Jo(k\r-r'\) (d = 2) (2.38a) 

Po 

i */ \ / iw u o(r) f +7T sin 6 'd9 , 

(ip (r)ip(r )} = / exp(ik\r — r \ cost)) 

Po Jo 2 

voir) sm(k\r — r'\) . , 

= -^^W JN 1 d = 3 . 2.386 

This equation is obviously valid only for \r — r'\ <C L, with L the typical size of the 
system. 

More generally, (2.37) with the notion that on average [<p]w(?,p) = fi mc (r,p) makes 
it natural to model the statistical properties of an eigenstate (pi close (on the scale of 
L) to some point r, by a superposition of a large number M of plane waves 

<p%( r ') = °v ex p(^ • ( r ' - r )/ h ) ( 2 - 39 ) 

ft=l,M 

where the p M are uniformly distributed with the probability 8(e n — H(r,p)) and, for 
time-reversal non-invariant systems, the Oj are complex random numbers such that 

1 Mr) 

M po 

Time-reversal invariance can be taken into account by having 2M plane waves (/x = 
±1, . . . , ±M) with = -pp, gl m = a*, and (a^a*,^,) = (v (r) /2p M)8 iV 5\ lx \\ lx >\. 

I shall for instance use this approach in section 6 to include the residual interactions 
into the fluctuations of Coulomb-blockade peak-spacings. In that case, the leading-order 
terms in an expansion in 1/g, where g = E Th /A the dimensionless conductance, can 
be derived straightforwardly from the model (2.39)-(2.40), supplemented only, following 
a kind of "minimum information hypothesis", by the assumption that the form a 
Gaussian vector. 

Sub-leading terms are however "aware" of the fmiteness of the system. One then 
needs to further modify the description of the wave-function fluctuations so that it also 
includes the typical scale L of the system under consideration. A simple way to do this 
is to use only plane waves fulfilling the quantization condition 

p^ = 2nhn^/L (2.41) 

with the set of integers n M = (ni, ■ ■ ■ , rid)- While doing so, one needs to give a width 
~ h/L to the Fermi surface, and to include in the modelling of the eigenstates all the 
plane wave with a kinetic energy in a band Se = aE^ around the Fermi energy, with 
a a constant of order one. For the M = ag basis vector in this shell, one can then 
use a random matrix description. To leading order in 1/M, the eigenstates fa then 
fulfil (2.39)-(2.40), as well as the Gaussian character of the coefficient a ifl , but some 
correlations of order 1/M are induced by the orthonormalisation constraints (Brody 
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et al. 1981). More importantly however, the width Sp ~ h/L of the Fermi surface 
modifies the wave-function correlations at distances of order L (for instance the two 
point correlation function decreases more rapidly than what is expressed by (2.38a)- 
(2.38&)). Time-reversal symmetry can here be taken into account using a real basis 
((cos(p M ■ r), sin(p M ■ r)) instead of (exp(±ip M • r)), and a GOE matrix. 

This model, which is defined both by the choice of the basis vectors and by the 
use of random matrices is what I will refer to as the random-plane-wave (RPW) model. 
It has to be borne in mind that this is a model, justified by physical considerations, 
but that should be in principle validated by comparing the fluctuations derived within 
this model to those obtained for the eigenfunctions of actual chaotic systems such as 
quantum billiards (see for instance in this respect McDonald & Kaufman 1988, Kudrolli 
et al. 1995, Pradhan & Sridhar 2002, Urbina & Richter 2003, Miller et al. 2005). 
In particular, relatively delicate questions concerning the normalization of the wave- 
functions may be important for some statistical quantities, which then requires 
further modifications of the random-plane-wave model described above (Urbina & 
Richter 2004, Urbina & Richter 2007, Tomsovic et al. 2008). 

2. 3. Screening of the Coulomb interaction in quantum dots 

For a degenerate electron gas, the "strength" of the Coulomb interaction 



between particle is usually expressed in terms of the gas parameter r s = r /oq, where 
a = f h 2 /m e e 2 is the (3c?) Bohr radius and r^ is the radius of a d- dimensional sphere 
containing on average one particle. Expressing, for instance for d = 2 or 3, the density 
of particle = (2irh)~ d g s J dp 0(e — p 2 /2m) as (g s = 2 is the spin degeneracy, and 
the Heaviside function) 



we see that, up to a constant of order one, 7q is essentially the inverse of the Fermi wave 
vector k F . As a consequence, r s is, again up to a constant of order one, proportional to 
the ratio (e 2 /r^) / '(ti 2 k F /2m e ) of the Coulomb energy between two electrons at typical 
inter-particle distance and the kinetic energy. The parameter r s therefore actually 
measures the relative strength of the Coulomb interaction. 

Even for small r s , U cou i(r — r') is long range, and for this reason large, in the sense 
that it cannot be taken into account by a low-order expansion. When physical properties 
are considered at an energy scale much smaller than the Fermi energy, it is however 
known (and well understood) that, for bulk systems, this interaction is renormalised 
because of screening into a much weaker effective interaction V sc (r— r'). Approximations 



Koui(r - r') 



r — r' 



(2.42) 




(2.43) 



(2.44) 
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for V sc (r— r') can be obtained using for instance the random phase approximation (RPA) 
(Fetter & Walecka 1971) giving in the zero- frequency low-momentum limit 

V sc (r) = J ^ V sc (q) exp[iq-r] , (2.45) 
2vre 2 

V sc (q) = -— (d = 2) , (2.46) 

\q\ + K (2) 



2 



V sc (q) = (d = 3) , (2.47) 

q -+- /t (3) 

with K{2) = (27re 2 )(g s z/Q 2 ' ) ) and k 2 ^ = (Aixe 2 )(g s h>^) the screening wave vectors + . 
One way to understand the screening mechanism is to view it in the spirit of the 
renormalization-group approach, where the effective Coulomb interaction that should 
be used for low-energy processes is produced by the integration of the "fast modes" 
(high energy degrees of freedom) of the electron gas (Shankar 1994). 

For finite systems, the situation is slightly more complicated because the 
renormalization process which transforms the bare interaction (2.42) into the screened 
one (2.46) or (2.47) also produces a mean field potential U m{ (r) which modifies the one 
particle part of the electrons Hamiltonian. Since both processes (screening and creation 
of the mean field potential) take place at the same time, their interplay is a priori not 
completely obvious. 

In the semiclassical limit, and more precisely whenever the screening length k," 1 is 
much smaller than the typical size L of the system, the common wisdom - that I shall 
follow here whenever necessary - is however simply to state that since the characteristic 
scales of variation of V^ c and of U m { are parametrically different (the former hT x is a 
quantum scale, when the latter L is classical), one could nevertheless use the same 
screened interaction as for the bulk, and furthermore assume that U m f(r) is correctly 
approximated by a Thomas-Fermi approximation. This latter amounts to minimise, 
with a fixed number of particles, the density functional 

J^tfM = T TF [n] + £ cxt [n] + £ coul [n] , (2.48) 

where 

e 2 f , , ,n(r)n(r') 
£couiN = 7T / drdr' V ; Y 
2 J \r — r'\ 

£ C xtM = J drU cxt (r)n(r) 

(£4xt( r ) is the external confining potential), and the kinetic energy term, originating 
from the Pauli exclusion principle, is given by 



TtfM = J drt TF (n(r)) 



+ The equations (2.46) and (2.47) are actually the expression of the screened interaction in the 
Thomas- Fermi approximation. To avoid confusion with the Thomas- Fermi approximation for the mean 
field potential, I will nevertheless refer to them, although slightly improperly, as the RPA-screened 
interaction. 
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t TF (n) = / dn' e(n') , 



(2.49) 



./o 

with e(n) the inverse of the function n(e) defined as 




(2.50) 



In particular e(n) = (h 2 /2m e )(4im/g s ) 2 for d = 2, et e(n) = (/l 2 /2m e )(67r 2 n/g s ) 2 / 3 for 
d = 3. 

The self-consistent equations obtained by minimising the Thomas-Fermi functional 
then read 



Note however, currently there is no general microscopic derivation of the above 
picture. More precisely, our confidence in having the Thomas-Fermi approximation 
as a correct starting point for the computation of U m f is due to the fact that this 
approximation can be derived in a quite general framework starting from a density 
functional description (in e.g. the local density approximation) and neglecting the 
effect of interferences (Ullmo et al. 2001, Ullmo et al. 2004). The "common wisdom" 
prescription given above therefore essentially amounts to trusting the density functional 
approach on the classical scale L (although it might be less reliable on the quantum 
scale Ap ; cf. for instance the discussion in Ullmo et al. (2004)), keeping the usual (bulk) 
form of the screened interaction on the quantum scale, and assuming that the two scales 
are not going to interfere in any significant way. That there is no microscopic derivation 
of this "common wisdom" prescription is presumably not too much of an issue as far 
as qualitative or statistical descriptions are concerned, but might become a limitation 
when accurate simulation tools are required to describe quantitatively the particular 
behaviour of a specific mesoscopic system. 

In this respect, one should note that there is a class of systems (namely billiards 
with weak disorder) for which it is possible to perform a renormalization procedure 
(Blanter et al. 1997, Aleiner et al. 2002) where the fast modes are integrated out so 
that only the interesting low-energy physics remains. It is then possible to see how 
both the mean field and the screened interaction emerge from this procedure. The 
generalisation to a more general case, for which £7 m f(r) is not well approximated by a 
constant, is however not completely straightforward, and is still an open problem. A 
relatively extensive discussion of this question is given in Appendix A. 

3. Orbital magnetism: general formalism 




(2.51) 



(2.52) 



As a first illustration of physical systems where the interrelation between interference 
effects and interactions plays a fundamental role, we shall consider in this section the 
orbital part - by opposition to the Zeeman part, associated with the coupling of the 
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magnetic field to the spin degrees of freedom - of the magnetic response at finite 
temperature k-^T = of mesoscopic objects in the ballistic or diffusive regime. 

Since the Bohr- van Leeuwen theorem (Leeuwen 1921), it is known that the 
magnetic response of a system of classical charged particles is exactly zero. This is 
a simple consequence of the fact that when writing the classical partition function 
Z c \ = J dpdrexp (—f3H c \{r,p)) with H c \(r,p) — (p — eA(r)) 2 /2m + U(r), the vector 
potential A(r), and thus any dependence in the magnetic field, can be eliminated by a 
change of the origin of p in the integral over momentum. The same holds true for the 
Weyl density of states since, up to an irrelevant multiplicative constant, it is derived 
from Z c \ by an inverse Laplace transform. As a consequence, whatever is measured has 
to be related to quantum effect \ , and, in the case of bounded fully coherent systems as 
is disordered systems, more specifically interference effects. 

In the early nineties, progresses in the design and probing of micron scale electronic 
systems, such as small metallic grains or quantum dots patterned in semiconductor 
hetero-structures, made is technically feasible to measure the magnetic response 
of coherent electronic structures. As the orbital magnetism is a particularly well 
adapted probe of interference effects in these coherent structures, this motivated a 
series of experimental work for (disordered/ diffusive) metallic grain (Levy et al. 1990, 
Chandrasekhar et al. 1991) as well as, slightly later, ballistic quantum dots (Levy 
et al. 1993, Mailly et al. 1993). 

All of these experiments were able to give convincing evidences - in particular 
the magnetic field scale - that the measured magnetic response was indeed related 
to quantum interferences. Moreover it was relatively soon realized that although the 
magnetic response of a single ring or dot can be dominated by terms for which the 
interactions are irrelevant, the dominant non- interacting contribution to the magnetic 
response was varying very rapidly with the size or chemical potential of the system. As 
a consequence, after averaging, the mean response of an ensemble of micro-structures 
is, most probably, driven by the contribution of the interactions. In other words, the 
magnetic response of ensembles of coherent electronic micro-structures is due to the 
interplay between interference effects and interactions. 

I stress however that "most probably" is the best that could be said here. Indeed, 
after the first series of experiments (Levy et al. 1990, Chandrasekhar et al. 1991, Levy 
et al. 1993, Mailly et al. 1993), which has sparked a host of theoretical works §, enough 
puzzles remain to indicate that a full understanding of the experimental data is lacking. 
I shall come back in more detail to this point at the end of section 4. It should be 
born in mind however that what follows should not be understood as the final "theory" 

| For instance, as discussed in Richter et al. (19966), the Landau diamagnetism can be understood as 
originating from quantum corrections to the Weyl term in the smooth density of states. 
§ See for instance (Bouchiat & Montambaux 1989, Ambegaokar & Eckern 19906, Ambegaokar & Eckcrn 
1990a, von Oppcn & Riedel 1991, Eckern 1991, Schmid 1991, Oh et al. 1991, Argaman et al. 1993, von 
Oppen & Riedel 1993, Gefen et al. 1994, von Oppen 1994, Ullmo et al. 1995, Montambaux 1996, Richter 
et al. 19966, Richter et al. 1996c, Ullmo et al. 1997, Ullmo et al. 1998, von Oppen et al. 2000). 
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of orbital magnetism in mesoscopic systems, but merely as the predictions that can be 
obtained for the equilibrium properties within a perturbative / Fermi-liquid description. 

In practice, I shall therefore consider a model of particles confined by some 
potential U(r), which is already assumed to contain the smooth part of the Coulomb 
interaction within some self-consistent scheme, interacting through the (residual) 
screened interaction (2.45), and subject either to a uniform magnetic field B, or, in 
the case of a ring, to a flux line $. Spin will be included only as a degeneracy factor 
(i.e. Zeeman coupling will not be considered). I shall furthermore limit the discussion 
to two-dimensional systems, but no drastically new effect is expected for d = 3. 

Within the grand canonical formalism, our goal will be to compute perturbatively 
in the interactions the field dependent part of the grand potential 

fi = -IlnZ G . c (3.1) 

with Zqc = Trexp(— (3(H — fiN)) the grand canonical partition function, and from 
there, for instance in the uniform field case, the magnetisation 

<*•> - -H (3 - 2 > 

(for completeness, equation (3.2), as well as its analog (4.26) for persistent current, is 
re-derived briefly in Appendix B) or the susceptibility 

x = a[-^b- ■ (3 - 3) 

(A is the area of the micro-structure). 

In the bulk, and more precisely when the cyclotron radius is larger than the 
coherence length or the thermal length Lt (to be defined more precisely below), 
the magnetic response is given by the (diamagnetic) Landau susceptibility Xl = 
— g s e 2 /24"7rm e . This latter originates from higher order in h corrections to the Weyl 
(smooth) density of states (Kubo 1964, Prado et al. 1994) (see also the discussion in 
section 3 of Richter et al. (1996&)), and I will use it below as the reference scale for the 
susceptibility. 



3.1. First order perturbations 

As mentioned in section 2.3, the effective strength of the screened interaction is related to 
the parameter r s characterising the density of the electron gas. In most experimentally 
relevant cases, r s is of order one and the high-density expansion is just a convenient 
way to order the various contributions, but some re-summation of series of higher order 
diagrams is necessary in order to get an accurate result. On the other hand, it is 
interesting and pedagogical to start with the genuine high-density asymptotics of small 
r s . Then, provided the momenta involved are of the order of the Fermi momentum 
Pf (which will be the case, except for the notable exception of periodic orbits, see 
discussion below) V sc (q) will be of order r s j vq and the diagrammatic development of the 
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thermodynamic potential is indeed a development in r s . In that case we are interested 
in the first order (or Hartree-Fock) correction in the screened interaction, which can be 
evaluated without drawing any Feynman diagram. 

3.1.1. First order perturbations Working in the grand-canonical ensemble at 
temperature k B T = f3 _1 , one can express the first correction to the thermodynamic 
potential as a direct (Hartree) plus an exchange (Fock) contribution in terms of 
the eigenfunctions ip u and eigenenergies e u of the non-interacting problem (Fetter & 
Walecka 1971), 

= 2 S fMSB 2 (<Pu<Pv\V\<p v <p v ) - g s ((p u <Pv\V\(p v (p u )] , (3.4) 

u,v 

with /„ = f(e v — /jl) = [1 + exp[/3(e„ — /i)]]^ 1 the Fermi occupation factor. For sake of 
clarity, the spin degeneracy factor g s = 2 is made explicit. 
Introducing 

n(r,r') = f v (r'\(p v )((p v \r) 

V 

= - ±J def(e-^[G R (r,r';e)-G A (r,r';e)], (3.5) 

(n(r) = n(r,r) is the local electron density in the non-interacting problem), one can 
re-express the direct and indirect contributions as 

U = 1 J drdr' n(r)V(r - r')n(r') (3.6a) 

F=\J drdr'n(r,r')V(r-r r )n(r',r) (3.6b) 

As discussed in Sec. 2.1, G R (r, r'; e) can be semiclassically approximated as a sum over 
classical trajectories travelling from r' to r at energy e. The advanced Green's function 
can be written in terms of the retarded one as 

G A (r,r';e) = [G R (r',r;e*)Y (3.7) 

and can therefore be interpreted as a sum running over all the trajectories going 
backward in time from r to r'. 

Using the equation (2.7) which relates the variation in energy of the action with 
the time of travel tj of the trajectory j, we understand the integral in (3.5) as the 
convolution between a function oscillating with a period 2irh/ tj and the Fermi function 
which varies from one to zero on a scale f3~ l = kBT. Introducing the characteristic time 
(or length) associated with temperature 

*r=^ = ^, (3.8) 

V F 7T 

(v F is the Fermi velocity) we see that the contribution of a trajectory j will be 
exponentially damped as soon as tj ^> t T . More precisely (see for instance the 
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appendix A of Richter et al. (19966)), if Gf is the contribution of the trajectory j 
to the sum (2.3), one has 

j de f(e - AGf(r, r') = (-^R( tj /t T )Gf (r, r')) , (3.9) 



with R(x) = f xj sinh(x). To proceed in the evaluation of the equations (3.6a)-(3.66), 
let us introduce the coordinates f = (r + r')/2 and 5r = (r — r'). Since the interaction 
between electrons is taken to be short ranged, one can assume the relevant 8r to be 
small and, using (dSj/dr') = pf, (dSj/dr) = —p), with p),Pj the initial and final 
momentum of the trajectory j, one can approximate 

iR(™ _i_ x~, In ™ _i_ x™' ln\ — nR 



Gf(f ± Sr/2, f ± 5r'/2) = Gf(f, f) exp 



.(3.10) 



The integral over 5r therefore yields the Fourier transform of V sc (r — r') and, neglecting 
the terms G A G A and G R G R which will eventually average to zero, one can semiclassically 
express the first order correction to the thermodynamic potential a sum over all pairs 
(k, I) of closed orbits 

H= {2^h I dr ^QkQi^s^ k -^)V sc ^p{-pl + p{-p})/2^) , (3.11a) 

T = J^Kfh j dr E^ C0S ^-^)^((^-^+pf-P/)/2^) , (3.116) 

where Qj = R(tj/tT)Dj/tj and tpj = (Sk/h— Qir/2). The field dependence of the above 
expression can then be obtained from (cf. 2.22) 

QQ. 

q£ = 2^/0o , (3.12) 

with &j the area enclosed by the orbit j. 

For a generic pair of trajectories (k,l) the term cos[(S' fc ± S t )/h] will be a highly 
oscillating function of the coordinate r. Performing the integration over position, the 
stationary phase condition reads {p{—p\) ± (pf — p\) = 0, and unless k and I are 
related by a symmetry, this will correspond to isolated points, each of which yields a 
contribution h 1 / 2 smaller than the original prefactor. On the other hand, special pairings 
where S k = Si will kill the oscillating phase. Such a condition is trivially satisfied when 
k = I but this also kills any field dependence. A second possibility is to pair a given 
trajectory with its time reversed. This is a nontrivial pairing since the resulting term 
is field dependent. Keeping only these contributions the direct and exchange terms can 
be written as 

/ HRitj/tr) \ 2 , _ , 2 _ / Air&jB\ y ( p) - P 



(3.13a) 
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where the sub-index (D) indicates the diagonal approximation, with sums running over 
individual trajectories j (and not pairs as in (3.11a)- (3.116)). A third possible pairing 
appears when we can match the actions of k and /, even if the trajectories are not the 
same or time-reversed of each other. Such situation arises in integrable systems with 
families of trajectories degenerate in action, and is discussed in details in Ullmo et al. 



3.1.2. The high density limit The sums in equations (3.13a)-(3.136) run over all closed 
(not necessarily periodic) trajectories [more precisely, we have twice the sum over all 
time-reversed pairs]. Therefore a priori, and in contrast to the non-interacting theory, 
periodic orbits do not play any particular role. It is interesting to note however that if 
the high density limit is to be taken seriously (i.e. r s — * 0), then again periodic orbits 
are singled out. Indeed, in this case V sc ([p'^—p^/fi) is of order except when 

p'j — pj = 0, that is when the orbit is periodic, in which case V^ c (0) = Vq 1 . Note that 
Pj + P j — implies that the trajectory j is self retracing, and thus has a zero enclosed 
area. As a consequence for the exchange term all contributions to the magnetic response 
are of order r s . 

Therefore in the high density limit, the integrand in (3.13a) is significantly larger 
in the neighbourhood of periodic orbits. For chaotic systems this will be compensated 
by the fact that these orbits are isolated. As a consequence the relative weight of 
their neighbourhood may depend on the particular system under consideration. For 
integrable systems however, for which periodic orbits come in families whose projection 
on the configuration space has a non-zero measure, the magnetic response induced by 
electron-electron interactions will be dominated by periodic orbits when r s « 1, and 
will reach a finite limit as r s — > 0. 

3.2. Correlations effects 

As discussed at the beginning of this chapter, realistic values of the parameter r s 
(appropriate for metals or GaAs/AlGaAs hetero-structures) force us to consider high- 
order effects in the diagrammatic expansion of the thermodynamic potential. Therefore, 
correlation effects are important, and need to be taken into account. 

A value of r s ~ 1 means that the range of the screened potential is of the order of 
the Fermi wavelength. In other words, the screened interaction has a local character, 
and can be written as 



where g s z/ is the total density of states (i.e. including the spin degeneracy factor g s ) 
and A is a constant of order one that will also serve to label the order of perturbation. 

The perturbative expansion of the thermodynamic potential can be represented 
in the usual way (cf. for instance section 15 of Abrikosov et al. (1963)) by Feynman 
diagrams, with straight lines standing for the finite temperature (Matsubara) Green's 



(1998). 




gs^O 



(3.14) 
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Figure 1. Direct terms of the Cooper series for the perturbative expansion of the 
thermodynamic potential f2. 

functions and the wavy lines for the interaction V . In the same way as in the theory 
of superconductivity, it turns out that it is essential to consider the full Cooper series 
which associated diagrams are represented in figure 1 (Abrikosov et al. 1963, Aslamazov 
& Larkin 1975). One way to see this is to realize that, since r s is actually not a small 
parameter, what is done here is more a semiclassical expansion (i.e. in powers of h) than 
one in the strength of the interaction. Therefore one should perform a simple counting 
of the powers of h of such contributions. The Cooper diagram of order k implies k 
interaction lines (each of which yields a factor Uq 1 ), n pairs of Green's functions, and 
(k + 1) summations of Matsubara frequencies (each of them associated with a factor 
fi^ 1 ). As far as powers of h are concerned, |G7| 2 ~ u /h (whatever the dimension). 
Therefore the only delicate point here is to realize that each temperature factor 
should be accounted for as an h, since in the mesoscopic regime considered here, the 
time tr = fr(3/ir introduced above should be of the order of some characteristic time 
t c of the system (say the time of flight), and thus ~ h/t c oc h (again as far as 
powers of fi are concerned). The RPA series can be seen, in the same way, to be of 
the same order in h, but the corresponding terms turn out to have negligible magnetic 
field dependence, and can therefore be omitted from the calculation of the magnetic 
response. Moreover one can convince oneself that all other diagrams would, at some 
given order k, have either a smaller number of Green's function or a larger number of 
frequency summations, and therefore are of higher order in h. 

Noting that, because the interaction (3.14) is local, the direct and exchange Cooper 
diagrams differ only by their sign and by a spin degeneracy factor, the magnetic response 
can be derived from the Cooper series contribution to the thermodynamic potential || 

fiC = * 28 E ~k E M r i • • • dr k Z(r u r k ; u m ) . . . S(r 2 , r i; u m ) 

k=l u) m <e F 

= ^3^ E Tr{ln[l + A S(r,r> m )]} , (3.15) 

^m<eF 

|| Note that the diagrammatic rules for ft differ slightly from the ones for correlation functions. There 
is in particular a factor 1/fc associated with each term of order k, thus the log. 



CONTENTS 



24 



where u m = 2nm/l3 are [bosonic] Matsubara frequencies, 

E(r, r'; cu m ) = — — V £(r, r'; e„)C?(r, r'; ^ m - e n ) (3.16) 

is the [free] particle-particle propagator, and the finite range ~ Af of the interaction 
introduces a cutoff on the summation over Matsubara frequencies at the corresponding 
energy scale e F . 

The trace over the space coordinates is a short way of expressing the expansion in 
all orders in AoS. The concept of particle-particle propagator, as well as the Cooper 
series contribution, comes from the Cooper pairs in the theory of superconductivity. 
The main difference in our case is that now the interaction is repulsive (thus the plus 
sign in the trace) and that we have lost translational invariance (therefore we cannot 
trade the operators for ordinary functions by going to the momentum representation). 



3.2.1. Semiclassical evaluation of the particle-particle propagator To proceed further 
with our semiclassical formalism, it is useful to write the finite-temperature Green's 
function between points r and r' for a fermionic Matsubara frequency (or rather, energy) 
e n = (2n + l)n/ P in terms of the retarded and advanced Green's functions as 

Q(r, r'; e n ) = Q(e n )G R (r, r'; e F +ie n ) + Q(-e n )G A (r, r'; e F +ie n ) , (3.17) 

with e F the Fermi energy. The retarded and advanced Green's functions are related 
through (3.7) and expressed, in a semiclassical approach, as expansions over all 
trajectories j joining r' and r at energy e (see (2.3)). The complex energy-arguments 
of (3.17) force us to perform some analytic continuation. However, if the Matsubara 
energies are much smaller than e F , one can expand the classical action and use (2.7) to 
obtain 



Gf(r, r'; e F +ie n ) = Gf(r, r'; e F ) x exp 



€ n tj 

"IT 



(3.18) 



Note that, as in (3.8) temperature introduces the time scale tx = h/3/ir which 
exponentially suppresses the contributions of long paths through the term e n tj/h = 
(2n + l)tj/t T . Therefore, only small Matsubara frequencies need to be considered, and 
the assumption used for the perturbative treatment of the energy is consistent. 

To compute the magnetic susceptibility at B = 0, the field dependent part of the 
semiclassical Green's function can also be treated perturbatively, and using (3.12) one 
can write 



Gf(r,r';e F +ie n ;B) = Gf(r,r';e F ;B = 0) x exp 



€ n tj 

'IT 



x exp 



B&j 



(3.19) 



where a^ is the effective area enclosed by the orbit (circulation of the vector potential 
between r and r') and 0o the flux quantum. The weak-field semiclassical approximation 
to (3.17) is then given by 



g(r,r';e n ,B)=6(e n 



1 E 

j:i — >t" 



, V-2inh 3 



x exp 



x exp 



i2n- 



B&i 



+ 
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(b) 



(d) 



Figure 2. Pairs of orbits contributing to Vl c (see (3.15) for a (non intcgrable) billiard. 
Top row : first order contributions. Bottom row : third order contributions (there are 
therefore three pairs of orbits connected at interaction points in both (c) and (d)). Left 
column : generic case. Right column : pairing of time reversed trajectories (diagonal 
contribution), for which the dynamical phases cancel. (Courtesy of Harold Baranger.) 



] -r ■ 



V-2inh 3 



-iS j/ /h+iTT( j ,/2 x 



cxp 



€ n tj' 



x exp 



-i2n- 



00 



(3.20) 



where trajectories j and f travel from r to r' in opposite directions, at energy e F , and 
in the absence of magnetic field. 

Note the usefulness of (3.20) goes beyond the problem of orbital magnetism 
discussed here, as it provides a calculational approach to any perturbative problem 
where the single-particle classical dynamics is known. 

The particle-particle propagator S(r,r';c<j m ) can now be evaluated semiclassically 
from (3.20) and (3.16). In general this involves a sum over all pairs of classical 
trajectories joining r' to r. An illustration is shown in figure 2. As in section 3.1 
however, most of these pairs yield highly oscillating contributions which average to zero 
when integrated over position, and one should only consider the non-oscillating terms 
which maintain a field dependence. One way to do this is, again, to pair time reversed 
trajectories, which implies that in the sum over the fermionic Matsubara frequencies in 
(3.16), only the e n such that e n (u m — e n ) < should be kept. This diagonal part 
of the particle-particle propagator can then be written as 

x exp 

(u!m — Cn)<0 

Summing over e n in the contribution of trajectory j, one gets 



\D-\ 2 



cxp 



(I* 



■A B ^3 



h 



)tj 



(3.21) 



e n (ui m -e n )<0 



UJ r , 



)U 



h 



= exp 



LO m tj 



R(2tj/t T ) 
2t 3 /t T 



cxp 



h 



,(3.22) 
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where the function R and the temperature time were introduced in the discussion of 
(3.9). The last factor (1 — exp[— (e F — u m )t/H]) originates from the upper bound e F of 
the Matsubara sum. If one assumes uj m <C e F , this factor removes from ^ D \r,r') all 
the contributions of trajectories of length smaller than A = X f /tt, thus preventing the 
particle-particle propagator from diverging as r — > r' . Replacing it by a hard cutoff at 
A one obtains 

g s H m e 2tj/t T 

Lj>A Q 

The semiclassical form for £( D )(r, r'; oo m ) shares with Hd and To (equations 
(3.13a)-(3.136)) the property of being a semiclassical expansion which does not oscillate 
rapidly (on the scale of \ F ) as a function of the coordinates, as would be the case for 
the Green's functions (3.20). 



iAti— 1 - 



exp[-— . 



(3.23) 



4. Orbital magnetism: Diffusive and ballistic systems 

4-1- Diffusive Systems 

The semiclassical approach described above does not rely on any assumption concerning 
the character of the underlying classical dynamics. It is therefore applicable to 
(integrable or chaotic) ballistic structures (Ullmo et al. 1998, von Oppen et al. 2000) 
as well as to diffusive systems (Montambaux 1996, Ullmo et al. 1997). Because 
diffusive motion is in some sense relatively simple, it is natural to consider first the 
orbital magnetism of interacting systems whose non-interacting classical dynamics is 
diffusive. More specifically, I will discuss the interaction contributions to the persistent 
current of metal rings and to the susceptibility of singly-connected two-dimensional 
diffusive systems. We shall see that in this case, the semiclassical approach recovers, 
in a very transparent and intuitive way, results previously obtained by quantum 
diagrammatic calculations (Ambegaokar & Eckern 19906, Eckern 1991, Aslamazov & 
Larkin 1975, Altshuler et al. 1983, Oh et al. 1991). Applied to diffusive dynamics, 
the semiclassical approach is indeed on the same level of approximation. Moreover, by 
making the connection with the classical dynamics, it provides a physically intuitive 
picture of the interplay between disorder and interactions. 

We assume here that the Fermi wavelength Xp is the shortest length scale, in 
particular \p < £ with £ the elastic mean free path, and that the magnetic field is 
classically weak, i.e., that the cyclotron radius at the Fermi energy is such that R c ^> £. 
Then the paths entering into (3.23) can be approximated by those of the system at zero 
field, but include the presence of the disorder potential. 

For diffusive systems it proves convenient to relate to the (classical) conditional 
probability P^(r,r';t\A) to propagate from r' to r in a time t and enclosing an area 
A since this probability satisfies a simple diffusion equation. For that purpose let us 
introduce an additional time and area integration 1 = f dt5{t — tj) f dA5(A — a^) in 
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(3.23) in order to make use of the sum rule (2.24), which, for a two dimensional kinetic 
plus potential Hamiltonian, and including the constraint on the area, reads 

E l ^^8(t-t j )6(A- f L j ) = 2nPZ(ry,t\A) ■ (4-1) 

j:r — >v 

One therefore obtains 

^ D \r,r';u m ) c f dA f dtP^r' \t\A)^Mxe W 

gs n> J Jt>A /v F zt /tT 
In the same way the n-th order (diagonal) contribution to the thermodynamic 
potential in (3.15) can then be expressed through the joint return probability 
P(t*i, r n , . . . , r*i; t n , . . . , h\A) to visit the n points r*j under the conditions that tj is 
the time elapsed during propagation from to T*j+i, and that the total enclosed area is 
A. For diffusive motion the probability is multiplicative, namely 

Jdr 1 ...dr n P(r 1 ,r n ,...,r 1 ;t n ,...,t 1 \A) = J drP(r,r;t tot \A) (4.3) 

with t tot = ^ti. Upon inserting the sum rule (4.1) and the relation (4.3) into (3.23), 
the contribution from the diagonal terms to f2 (see (3.15)) yields 

^ (D) = E^ D) 

n 

= ]l l ^ l dtC ° th {j~) K i t ) A i r ^ B ) ■ ( 44 ) 

The factor coth(t/t r ) (with t T defined by (3.8)) arises from the tu-sum in (3.15) which 
here can be performed explicitly. In (4.4) the functions K and A are defined as 



exp[-^] .(4.2) 



dUR(2U/t T ) 



5(t-t tot ) } ,(4.5) 



K(t) =Y. K ^) ; = { /n 

n \ i=l 

A(r,t;B)= J dAcos (^^j P(r,r-,t\A) . (4.6) 

K(t) accounts for temperature effects while A(r, t; B) contains the field dependence and 
the classical return probability. 

4-1.1. Renormalization scheme for Diffusive Systems This semiclassical approach 
allows us further to obtain the renormalization of the coupling constant (Altshuler 
et al. 1983, Altshuler & Aronov 1985, Eckern 1991, Ullmo et al. 1997) for diffusive 
systems by resuming the higher-order diagrams of the Cooper series. To this end let us 
first introduce the Laplace transform of Ki(t), 

Hp) = ^t ptT+2 \ 2n+iy (4.7) 



where 



f3e F k F L T 
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The full kernel K(t) is given by the inverse Laplace transform 

K(t) = — / dpe^Hl + fip)}. (4.9) 



2ni 

To evaluate the above integral, let us define 

g(p) = l + f(p), (4-10) 
and furthermore denote the singularities of g(p) by 

v. - (4 . n) 

It 

with n = 0, . . . , n F . Let p n be the corresponding zeros (p n is assumed to lie between p n 
and Pn+i)- On the real axis, g is a real function which is negative within each interval 
\pn,Pn] (with the notation p np = — oo) and positive elsewhere. As a consequence, lng(p) 
is analytic everywhere in the complex plane except for the branch cut \p n ,Pn]- The phase 
jump across the branch cuts is 2n, since the imaginary part of g(p) is positive above 
and negative below the real axis. Deforming the contour of integration as sketched in 
figure 3 one therefore finds 

K(t) = lim / ^- (ln[£(p - %e)\ - ln[£(p + %e)\) e pt (4.12) 

= X / ( 4 - 13 ) 

n=0 "'P" 
ra=0 

For n <^ n F one has 5„ = trG°n — Pn) ^ 1 an d thus to first order in 5 n : 

1 + - E ~~ T~ — " — ^ = • ( 4 - 15 ) 
gs^^'-n g s S n 

The above condition gives 

(4.16) 



g s /Ao + ^(n F + 1) - <S>(2n + 1) 
with \I> the digamma function. 

In the high temperature regime t T < t, all the n's actually contributing to the sum 
(4.14) are such that the denominator in (4.16) is dominated by $(rif + 1) ~ ln(np). 
One obtains in this case 

n F 

K(t) = - E ePnt i 1 - e~ 5 " t/tT ] (4.17) 



t 

n=0 



np 



In the low temperature regime tr/t ^> 1, the typical n contributing to (4.14) is 
no = t T /At (that is still assumed much smaller than rip). Because of the slow variation 



CONTENTS 



29 




Figure 3. Top: graph of the function g(p) (for n/ = 200). Bottom: integration path 
for the inverse Laplace transform (4.13) in the complex p plan. 



of the logarithm, one can in this case replace n by n in (4.16). [A formal derivation 
would consist in replacing the sum (4.14) by an integral, change to the variable ln(n), 
and use the stationary phase approximation.] This gives 

4 4 

6n ~ ln{n f ) - ln(2n ) = \n{2k F v F t) ' (4 ' 19) 
and in the same way as above 



n F 



K(t) ~ — — r Y e Pnt . (4.20) 

v ; tr \n(2k F v F t) v ; 

v ' n=0 
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Noting that 



4A np 

K ^) = vrY. ePnt ( 4 - 21 ) 



one identifies 



with 



min(2v F t, L T /4) . (4.23) 



This relation is valid when \nk F L* » 1. It is certainly satisfied when ln/c F / 3> 1 which 
may be regarded as a semiclassical approximation in the diffusive regime. The equation 
(4.22) shows that the higher-order terms in K(t) merely lead to a renormalization of 
the first-order contribution Ki(t): The coupling constant A = 1, entering into Ki(t), is 
renormalised to g s / ln(/c F L*). 

In first order one has from (4.5) 

K l( t) = (4.24) 
so that K reduces to 

Equations (4.4)-(4.6) together with (4.25) may serve as a general and convenient starting 
point to compute the orbital response of disordered mesoscopic systems. The specific 
character of the geometry enters into QS D ^ solely through the Fourier transform (4.6) 
of the return probability P(r, r,t\A). In the following we shall apply this approach to 
compute the magnetic response of two important types of mesoscopic structures. 

4-1-2. Persistent current of disordered rings While the magnetic response of a singly- 
connected system is usually described in terms of its susceptibility, the response of a 
ring-type structure threaded by a flux is related to the persistent current 

<«■*> 

(see Appendix Appendix B for a derivation of this equation). To make contact with 
previous approaches (Ambegaokar & Eckern 19906, Montambaux 1996) let us start with 
the computation of the contribution of the first-order interaction term, Q[ D ^ , to / which 
will be denoted by I±. 

Consider a (thin) disordered ring of width b, cross section o and circumference L. 
For L ^> l,b the motion of particles around the ring effectively follows a law for one- 
dimensional diffusion. The total area enclosed by a path is given in terms of the number 
m of windings around the ring. One thus has 
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where D = v F l/d is the diffusion constant in d dimensions and v F the Fermi velocity. 
Note that due to the disorder average, the classical return probability does not depend 
on r. 

Combining the expression (4.24) for Ki(t) with the coth function in (4.4) one 
obtains 

^ = ^?-(T)- m (4 ' 28) 



with 



9m(T) = (" **<fM «p l-MjV^)] ( 

Jo < VAnDt 

After taking the flux derivative according to (4.26), one recovers the first-order 
interaction contribution to the persistent current (Ullmo et al. 1997), 



& * m ~oo V A> 7 
This first-order result was first obtained by purely diagrammatic techniques by 
Ambegaokar & Eckern (1990&) and semiclassically by Montambaux (1996). 

However, higher-order terms are essential for an appropriate computation of the 
interaction contribution. As shown in the preceding subsection, these higher-order 
terms merely lead to a renormalization of the coupling constant according to (4.22). 
Thus the persistent current from the entire interaction contribution is reduced to (Ullmo 
et al. 1997) 

For diffusive rings the length scale v F t, entering in (4.23) for L*, is given by L m = 
v F (mL) 2 /4:D, the average length of a trajectory diffusing m times around the ring. 
Hence one gets at low temperature (L T ^> L m ) a (renormalised) prefactor ~ 2/ ln(/c F L m ) 
for /. At high temperature, Lt <C L m , the prefactor includes 2/ln(/c F Ly/4). These two 
limits agree with quantum results obtained diagrammatically by Eckern (1991). The 
functional form of the temperature dependence (exponential T-damping (Ambegaokar 
& Eckern 19906)) is in line with experiments (Levy et al. 1990, Chandrasekhar 
et al. 1991, Mohanty et al. 1996). However, the amplitude of the full persistent current 
with renormalised coupling constant is a factor ~ 3—5 to small compared to experiments. 

4-1.3. Susceptibility of two-dimensional diffusive systems In ring geometries the 
exponential temperature dependence of / is related to the existence of a minimal length, 
the circumference of the ring, for the shortest flux-enclosing paths. In singly connected 
systems the geometry does not constrain returning paths to have a minimal length, and 
therefore one expects a different temperature dependence of the magnetic response. 

Consider a two-dimensional singly connected quantum dot with diffusive dynamics. 
If one makes use of the general renormalization property of diffusive systems, expressed 
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by (4.25), the diagonal part of the thermodynamic potential (see (4.4)), including the 
entire Cooper series, reads 

The parameter L* appearing in (4.22) has been here replaced by v F t since the factor R 2 
ensures that the main contribution to the integral comes from t < tj 1 - 111 the above time 
integral the elastic scattering time r e i — l/v F enters as a lower bound. This cutoff must 
be introduced since for backscattered paths with times shorter than T e \ the diffusion 
approximation (4.33) no longer holds. Short paths with t < r e i, which may arise 
from higher-order interaction events, contribute to the clean bulk magnetic response 
(Aslamazov & Larkin 1975, von Oppen et al. 2000). This latter term is, however, 
negligible compared to the disorder induced interaction contributions considered here. 

To evaluate the integral for A, the conditional return probability in two dimensions 
is conveniently represented in terms of the Fourier transform (Argaman et al. 1993) 

P(r,r,t\A) = — / dk \k\ — tttttttt • ( 4 -33) 
v ' ' 1 ; 8tt 2 J 11 smh(\k\Dt) v ; 

Introducing the magnetic time 

* = *ss = £b • ^ 

(ts is related to the square of the magnetic length L 2 B which can be regarded as the 
area enclosing one flux quantum (assuming diffusive dynamics)), one obtains 

The function R occurring in (4.35) has a different origin than in (4.32). 

The magnetic susceptibility (3.3) is obtained after including the expression (4.35) 
in (4.32) and taking the second derivative with respect to the magnetic field. One finds, 
after normalization to the Landau susceptibility of non-interacting particles in a clean 

system, 

f. = J2 {krl) r * # (L) g. (1) . (4.36) 
vr V J Tel tln{k F v F t) \t T J \t B J 

Here, D = v F l/2, and R" denotes the second derivative of R. 

The above equation (4.36) holds true only as long as the effective upper cutoff 
time t* = min(tT , t b) , introduced through the i?-functions, remains smaller than the 
Thouless time t c = L 2 /D with L being the system size. For times larger than t c the 
two-dimensional diffusion approximation is no longer valid since the dynamics begins to 
behave ergodically. 

Assuming r e ; <C t* < t c , the integral in (4.36) can be approximately evaluated by 
replacing the upper bound by t*, and by replacing R{t/tr) and R"{t/tB) by R(0) = 1 
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and R"(0) = —1/3, respectively. The averaged magnetic susceptibility of a diffusive 
two-dimensional quantum system reads under these approximations (Ullmo et al. 1997) 



The magnetic response of diffusive systems is paramagnetic and enhanced by a factor 
k F l compared to the clean Landau susceptibility xl- 

Contrary to the exponential temperature dependence of the ring geometry discussed 
in the previous section one finds a log- log temperature dependence for tx < ts as well 
as a log-log B dependence for tr > ts- The log-log form of the result is produced by 
the 1/tlnt dependence of the integral in (4.36). It results from the wide distribution of 
path-lengths in the system — there are flux-enclosing paths with lengths ranging from 
about v F r e i up to v F t*. 

The expression (4.37) agrees with results from Aslamazov and Larkin (1975), 
Altshuler, Aronov and Zyuzin (1983), and Oh, Zyuzin and Serota (1991) which are 
obtained with quantum diagrammatic perturbation theory. The equivalence between 
the semiclassical and diagrammatic approaches, demonstrated here as well as the 
preceding section for diffusive rings may be traced back to the fact that the "quantum" 
diagrammatic perturbation theory relies on the use of the small parameter l/k F £, and 
can therefore also be viewed as a semiclassical approximation. 

4-2. Ballistic Systems 

We turn now to the magnetic response of ballistic mesoscopic objects. This problem 
has received considerably less attention than the diffusive regime, partly because of 
the experimental difficulties involved, but also presumably because from the theoretical 
point of view, ballistic systems cannot be addressed with the more traditional approaches 
of solid state physics. Indeed, for diffusive systems, the main virtue of the "quantum- 
chaos" based approach presented in this review was to provide an alternative, maybe 
more intuitive, way to get known results. In contrast, the ballistic regime represents 
one example for which it is actually the only way to get understanding of the physics 
involved. 

The ballistic character of the underlying classical dynamics brings about some 
important differences, especially since one no longer talks about probabilistic concepts 
(like the return probability), but one needs to input the information about specific 
trajectories. Depending on the geometry of the confining potential, one can have chaotic 
or integrable motion, and the structure of the semiclassical expansions differs in some 
respects in these two cases. In particular, the existence of families of trajectories for 
integrable dynamics will translate into larger magnetic response. 

4-2.1. Renormalization scheme for ballistic systems As for diffusive systems, (3.15) and 
(3.23) form the basic expressions from which the magnetic response can be computed. 
We just saw that in this latter case, they allow quite straightforwardly to express the 




(4.37) 
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magnetic susceptibility of dots or persistent current of rings in closed form. For ballistic 
systems on the other hand, even when the classical dynamics is simple enough to yield 
E( D '(r, r'\ u) m ) explicitly, taking the logarithm of this operator as implied by (3.15) will 
require proceeding numerically. It should be born in mind however that this numerical 
calculation is significantly simpler than the one we would have faced had we decided 
to start directly with a numerical approach. This can be understood for instance by 
considering that the original operator £(r, r'\ uj m ), because of its quantal nature, has a 
scale of variation determined by X F . Discretizing this operator on a grid would therefore 
require using a mesh containing at least a few points per X F , which would make any 
computation rapidly intractable as the size of the problem increases. On the other 
hand, within the semiclassical approach, the operator one has to deal with is £( D )(r, r') 
which, up to one exception to which I shall return below, varies only on a classical scale. 
£( D )(r*, r') is therefore what one may call a "classical operator". It permits the use of a 
grid mesh whose scale is fixed by the classical dynamics within the system, and therefore 
much larger than \ F . 

As mentioned above Y^ D \r, r') is not yet completely classical. Indeed, there still is, 
in (3.23), the scale A = X F /ir which specifies that trajectories shorter than A should 
be excluded from the sum over trajectories joining r' to r. This last quantum scale 
can be removed using a simple renormalization scheme where integration over short 
trajectories yields a decreased effective coupling constant. To that end consider a new 
cutoff A larger than A but much smaller than any other characteristic lengths (a, L?, 
or a/0o/-B)- For each path j joining r' to r with Lj > A, let T,j(r,r') denotes its 
contribution to T^ D \r,r') and define 

£,(r,r') = Ej(ry) - X j dr, E,(r, ri)E(n, r') 

+ X 2 J dr x dr 2 E,(r, r 2 )E(r 2 , r 1 )E(r 1 , r') + ■ ■ ■ . (4.38) 

where the r; integration is over A < |r*i_i — r*i| < A (with r = r'). E(?*i, r') is defined 
by (3.23) but with the sum restricted to "short" trajectories with lengths in the range 
[Ao,A]; ^j( r > r i) i s obtained from Ej(r, r') by continuously deforming trajectory j. To 
avoid the awkward In in (3.15), let us introduce 

1 



-YTt 



l + AoSW^r';^ 



(4.39) 



from which can be derived through 

nW(Ao)= / ^r(Ai). (4.40) 

The replacement E by E in T amounts to a reordering of the perturbation expansion of 
T in which short paths are gathered into lower-order terms. Moreover, if Lj 3> A, small 
variations in the spatial arguments do not modify noticeably the characteristics of Ej. 
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Approximating £,-(r, t*i) by £,-(r, r') in (4.38) and using £( r i; r ') — l/47r|r! — r'| 2 valid 
for short paths, one obtains 

A S,(r,r') ~ -^fiff - ^ A(A)E,(r,r') (4.41) 
1 + A J driE(ri,r') 

where the running coupling constant is defined by 

Therefore, these successive steps amount to a change of both the coupling constant and 
the cutoff (since now trajectories shorter than A must be excluded) without changing 
T; that is, 

r(A ,A ) = r(A,A(A)) . (4.43) 

Through (4.40), this renormalization scheme can be applied to Q( D \ and so to the 
average susceptibility. 

In this way, we have eliminated the last "quantum scale" Ao from the definition 
of £( D ); A can be made much larger than A F while remaining smaller than all classical 
lengths. This will serve two purposes. From a quantitative point of view, it means that 
discretizing the operator £( D ) can be done on a relatively coarse grid, and that therefore 
operations such as taking the trace logarithm of 1 + S^ D ^ that may be necessary to get 
the magnetic response can be done numerically much more easily. 

Furthermore, it is qualitatively reasonable that the perturbation series of 
becomes convergent when A is of order of the typical size L of the system, since by this 
point the spread in length scales causing the divergence has been eliminated. Of course 
the renormalization scheme above assumes, A much smaller than all the classical scales 
of the system, and using A ~ L is beyond the range for which a reliable quantitative 
answer can be obtained. At a qualitative level however, this implies that higher-order 
terms in the diagonal contribution mainly renormalise the coupling constant Ao into 

= i + (V gW) ■ (4 ' 44) 

In the deep semiclassical limit \n(k F L) ^> 1, the original coupling constant Ao drops 
from the final results and is replaced by g s / \n(k F L). 

For integrable systems in which non-diagonal channels exist, simple inspection 
shows that these latter are not renormalised. Indeed the corresponding higher- 
order terms in the perturbation expansion are highly oscillatory. The non-diagonal 
contributions are therefore of order Ao rather than g s / ln(k F L), and will as a consequence 
dominate the magnetic response in the deep semiclassical regime. 



4-2.2. Ballistic squares As an illustration, let us consider square quantum dots. The 
choice of squares is motivated first because this is the geometry used experimentally 
(Levy et al. 1993), and secondly because squares are particularly amenable to a 



CONTENTS 



36 



r 



Figure 4. Typical pairs of real-space trajectories that contribute to the average 
susceptibility to first order in the interaction in the diagonal channel (left) and the 
non-diagonal channel (right). This figure is taken from Ullmo et al. (1998). 

semiclassical treatment, since it is very easy to enumerate all the classical trajectories \. 
To be more specific, I will consider a square billiard model, of size L, with free motion 
(i.e. just a kinetic energy term) inside the the billiard and Dirichlet boundary conditions 
on its border. An illustration of the shortest (with non-zero area) closed orbit and of 
the shortest periodic orbit, for the classical version of the billiard, is given in figure 4. 

Let us discuss first the diagonal contribution Q( D ^ (3.15) to the thermodynamic 
potential. Because of the simplicity of the geometry, the explicit expression of 
Hj(r,r';u) for any orbit j joining r' to r can be obtained quite straightforwardly 
(though the resulting expressions might be a bit cumbersome and will thus not be 
given). At a given temperature T, the operator £( D )(r, r'; uo) is then constructed by 
summing all such contributions for orbits of lengths shorter than the thermal length 
Lt (see (3.8)). A numerical computation of (and therefore, after derivation with 
respect to the magnetic field, of the magnetic susceptibility) can then be obtain by 
representing S^ D ^(r, r'\ to) on a grid, then going to a diagonal representation, and in 
this diagonal representation taking the log and performing the trace. The temperature 
dependence of the resulting contribution to the susceptibility, for a number of particles 
corresponding to k F L = 50, is shown as a dotted line in figure 5. 

In this way, one can furthermore check that the renormalization process described 
above actually works for the square geometry. For instance figure 6 shows the value of 
the magnetic susceptibility at zero field, and for a temperature such that Lt = AL, as 
a function of the cutoff A. As long as A <C L, the diagonal part of the susceptibility is 
clearly independent of the cutoff if the renormalised interaction A(A) given by (4.42) is 
used. 

| As a consequence of the two previous facts, orbital magnetism in ballistic squares has been thoroughly 
studied within non-interacting models (von Oppen 1994, Ullmo et al. 1995, Richter et al. 1996&, Gefen 
et al. 1994, Richter et al. 1996c, Richter et al. 1996a). 
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Figure 5. Temperature dependence of the zero-field susceptibility (solid line) for an 
ensemble of squares at k F L = 50. The contribution of the non-diagonal channel (dashed, 
family (11) and repetitions) exceeds that of the diagonal Cooper channel (dotted) at 
low temperatures. Temperatures are expressed in units of k^T a = hv F /2nL. Inset: 
same, but with a different scale for the vertical axis. This figure is taken from Ullmo 
ct al. (1998). 



Let us finally compare the diagonal and non-diagonal contributions. This latter 
is built from pairs of orbits with the same action (i.e. for a billiard, same length) but 
different geometry, overlapping at a given point. Those may occur quite generically for 
integrable systems for which families of periodic orbits exist and two members of the 
same family may intersect in configuration space. An illustration of such an intersection 
is given in the right panel of figure 4 for the family (labelled (1, 1)) of shortest periodic 
orbits (with non-zero enclosed area) of this system. In the square billiard, for a given 
closed orbit j of length Lj and enclosing the action in the absence of magnetic 

field is given by Sj/h = k F Lj, its derivative with respect to the magnetic field by 
dSj/dB = 2na,j/(j) , and the stability determinant (cf. (2.5)) by Dj = m/ ^Jhk F Lj. The 
contribution °f the pairs of orbits from the family (1,1) can for instance be derived 
quite straightforwardly from equations (3.11a)-(3.11a) noting that that L n = 2\[2L and 
a n = 2x (L — x ) (x is the abscissa of the intersection of the orbit with the bottom 
border of the billiard, and can be use to label an individual orbit within the family 
(1, 1)). One obtains in this way (Ullmo et al. 1998) 

(X(T,i) diag ) Sk F L d 2 C 2 ((p) 2 f L n 



, , R I J . (4.45) 
\Xl\ 2(v^tt) 3 \L T ) K J 

The temperature dependence is governed by the function R{x) — x/ sinh(x) and the field 
dependence by C((f) = (2(p)~ 1 ^ 2 [cos(7T(/9)C( v /7t^) + sm(7^^p)S(^/K(p)] , with ip = BL 2 /(f) , 
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(B=0, k F L=50, L T =4L) 
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Figure 6. Diagonal part of the interacting contribution to the magnetic susceptibility 
as a function of the cutoff A used to define the particle-particle propagator S. The data 
correspond to zero magnetic field, k F L = 50, and a temperature T such that Lt = 4Z>. 
Solid line: computation for which the coupling constant A varies as a function of A 
following (4.42). Dash line: A is kept fixed (= Ao). The vertical (dash) line corresponds 
to the limit of the range for which the condition A <C L is fulfilled 



and C and S Fresnel functions. 

As seen in figure 5, because the non-diagonal contribution is not renormalised by 
higher-order terms, it dominates the magnetic response as soon as the temperature is 
low enough not to suppress it. 

4-3. Discussion 

After this overview of some of the results that can be derived for the magnetic 
response of mesoscopic rings or dots within a Fermi-liquid, at equilibrium, description, 
it is worthwhile to come back to how well existing experiments can be understood 
within this framework. Since the typical response for a single system (Mailly 
et al. 1993, Chandrasekhar et al. 1991) may in some case be dominated by the non- 
interacting contribution, which is not addressed here, I will however limit this discussion 
to experiments made on ensembles of micro-structures. 

Ballistic and diffusive systems differ with respect to the comparison between 
experimental measurements and theoretical predictions, but for both cases there is 
clearly not a complete adequacy. 

In the case of ballistic systems (Levy et al. 1993), the amplitude of the magnetic 
response in the low temperature range, as well as scale of the field dependence, are well 
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in line with the predictions (von Oppen 1994, Ullmo et al. 1995, Richter et al. 1996&). 
However, the temperature dependence, and in particular the fact that a magnetic 
response is observed even when the thermal length Lt is smaller than the size of the 
system seems very difficult to interpret within the Fermi liquid framework used here 
(Ullmo et al. 1995, Ullmo et al. 1998). There is, up to now, no real suggestion to explain 
how a significant magnetic response, with a field scale unambiguously associated with 
interference effects, could survive at such high temperatures. 

For diffusive rings (Levy et al. 1990), the temperature dependence seems to be 
less of an issue than in the ballistic case. However the magnitude of the average 
persistent current appears to be quite a bit larger (a factor three to five) than what 
is expected theoretically (Ambegaokar & Eckern 1990a). It has been suggested that 
this large magnetic response might be related to non-equilibrium effects, such as the 
coupling to a phonon or photon bath (Mohanty 1999, Kravtsov & Altshuler 2000, Entin- 
Wohlman et al. 2003). In later experiments (Jariwala et al. 2001) (see also Reulet 
et al. (1995) and Deblock et al. (2002) in this context), a further, and presumably 
more dramatic, discrepancy is that the sign of the magnetic response - which has not 
been determined in Levy et al. (1990) - was furthermore not the one expected for a 
repulsive interaction. This has motivated consideration of whether a BCS interaction, 
weak enough to be compatible with the absence of observed superconducting state for 
the considered materials, could be responsible for both the change of sign (which would 
be an immediate consequence of having an attractive interaction) but also the increase in 
the magnitude of the observed signal (Schechter et al. 2003, Eckern et al. 2004, Schechter 
et al. 2004). It is not clear however that such a range of attractive interaction actually 
exists in practice. 

In the end, it is remarkable that on the one hand a good part of the experimental 
findings has a simple and natural interpretation within the "confined electron gas with 
weak interaction" picture developed here, when on the other hand no consensus has 
emerged yet concerning the interpretation other significant experimental findings. One 
of the limiting factors in this respect is presumably the imbalance between the large 
number of theoretical works and the much sparser character of experimental studies on 
these questions. In particular, only one group has measured the magnetic response of an 
ensemble of micro-structures in the ballistic regime (Levy et al. 1993). In that case the 
low temperature data agree very well with the theoretical description when there is little 
hope to interpret what is observed by the high temperature regime without introducing 
some new physical mechanisms. However, the lack of experimental indications on what 
these other mechanisms could be is making difficult any theoretical progress. In the 
same way, for diffusive systems, a complete understanding would be made easier if one 
knew how much the sign of the magnetisation observed in Jariwala et al. (2001) is 
material dependant. 

Obviously, it would be very non-trivial to improve on the existing experiments, 
in particular because it is necessary to deal properly with very weak magnetic fields. 
However, there are enough indications that some deep physics might be involved to 
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motivate further experimental work in this field. 
5. Mesoscopic Kondo effect 

Turning now to another context where the interrelation between interferences and 
interactions plays an important role, I shall wander briefly away from Fermi liquids, and 
consider some (limited) aspects of the problem of a Kondo impurity within a mesoscopic 
conductor. 

5.1. A quick background 

The term "Kondo effect" refers to the physics of an impurity with some internal degree 
of freedom, interacting with a gas of otherwise non interacting electrons. It represents 
one of the simplest models in condensed matter physics for which correlations play a 
central role (Hewson 1993). 

In its simplest version, the s-d model, the impurity is just treated as a spin one half 
interacting locally with the electron gas. The corresponding Hamiltonian then reads 



where & aa creates a particle with energy e a , spin a and wave-function v?a( r )> and the 
interaction with the impurity is expressed as 



with J the coupling strength, S = (S x , S y , S z ) a spin operator (h 1 Si is half of the Pauli 
matrix a t ), s(0) = f (0)<w# ff (0) the spin density of the electron gas at the impurity 
position r = 0, and tft(0) = £ a y? a (0)ct. 



Originally, physical realizations of the Kondo Hamiltonian corresponded to actual 
impurities (e.g. Fe) in a bulk metal (e.g. Cu). The wave-functions tp a could then be 
taken as plane waves, and one could assume a constant spacing A between the e a , so 
that the electron gas could be characterized by only two quantities: the local density 
of states vq = (AA)^ 1 (A is the volume of the sample), and the bandwidth D of the 
spectrum. 

What gave (and still gives) to the Kondo problem its particular place in condensed 
matter physics is that it is the simplest problem for which the physics is dominated by 
renormalization effects. Indeed, assuming the dimensionless constant J u 1 , it can 
be shown using one-loop renormalization group analysis (Fowler & Zawadowzki 1971), 
or equivalent earlier approaches such as Anderson poor's man scaling (Anderson 1970) 
or Abrikosov parquet diagrams re-summation (Abrikosov 1965), that the low energy 
physics remains unchanged by the simultaneous modification of both Jo and Do into 
new values J e g and D cff provided they are related by 




(5.1) 



#int = T§ S ■ S(0) 



(5.2) 



Jeff(Dcflf) 



1 - J z/ hi(D /D eff ) ' 



(5.3) 
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The renormalization procedure should naturally be stopped when D e g- becomes of the 
order of the temperature T of the system. Equation (5.3) defines an energy scale, the 
Kondo temperature 

T A: = Doex P (-l/(Jo^o)), (5-4) 

which specifies the crossover between the weakly and strongly interacting regime. For 
T 3> T K) the impurity is effectively weakly coupled to the electron gas, and the 
properties of the system can be computed within a perturbative approach provided 
the renormalised interaction J e s(T) is used. The regime T <C Tk is characterized by 
an effectively very strong interaction (in spite of the bare coupling value J uq being 
small), in such a way that the spin of the impurity is almost completely screened by the 
electron gas. Perturbative renormalization analysis (and thus (5.3) itself) can obviously 
not be applied in this regime, but (for bulk system) a rather complete description 
has been obtained by a variety of approach, including numerical renormalization group 
(Wilson 1975), Bethe anzatz techniques (Andrei 1980, Wiegmann 1980), and in the very 
low temperature regime Nozieres Fermi liquid description (Nozieres 1974). 

One important consequence of the scaling law (5.3) is that physical quantities can 
be described by universal functions, which can be understood by a simple counting 
of the number of parameters defining the s-d model in the bulk. I will illustrate this 
discussion with a particular physical quantity namely the local susceptibility 

rf> 

Xioc = / dr(S z {r)SM) , (5-5) 
Jo 

which is the variation of the impurity spin magnetisation to a field applied only to 
the impurity. The electron gas is characterized by its local density of states v$ and 
its bandwidth D , and the impurity by the coupling constant J . Therefore, for a 
given temperature T, xioc, as any physical quantity, can depend only on these four 
parameters. Furthermore, only two dimensionless parameters can be constructed from 
them, the ratio T/D and the product JqVq- However, because of the scaling law (5.3), 
we see that these two parameters turn out to be eventually redundant. A dimensionless 
quantity can therefore be expressed as a function of a single parameter, which can be 
chosen to be T /Tk so that we have for instance 

Txioc = f x (T/T K ) , (5.6) 

with f x (x) a universal function which has been computed by Wilson (1975) using 
his numerical renormalization group approach. Note finally that within the one-loop 
approximation (5.3)-(5.4), 

T/T K = exp(J cS (T)v ) (5.7) 

so that another way to express the universal character of physical quantities is to say 
that they depend only on the dimensionless quantity J c q(T)u q . 



CONTENTS 



42 



5.2. Mesoscopic fluctuations 

We see that the universal character of the Kondo physics is a direct consequence of 
the fact that the local density of states is flat and featureless, and can therefore be 
characterized by two parameters, u and D . There are many situations however for 
which the variation in energy (and position actually) of the local density of states 
z/ ioc( r ;e) might be significantly more complex, and it is natural to ask in which way 
this would modify the description given above in the bulk (flat band) case. 

The origin of fluctuations in the density of states can be of different nature, 
but they are always in the end associated to interference or finite size effects. One 
possibility is the presence of disorder in a bulk material (Dobrosavljevic et al. 1992, 
Kettemann 2004, Kettemann & Mucciolo 2006, Kettemann & Mucciolo 2007, Zhuravlev 
et al. 2007), in either the metallic or localised regime. Another is the proximity of some 
boundary in the host material (Ujsaghy et al. 2001), for instance in the case of a narrow 
point contact (Zarand & Udvardi 1996) or for thin film (Crepieux & Lacroix 2000). 
Finally the class of systems where a Kondo impurity is placed within a fully coherent, 
finite size electron has been realized for instance in the context of "quantum 

corrals" (Fiete et al. 2001), has been also considered (Thimm et al. 1999, Affleck 
& Simon 2001, Cornaglia & Balseiro 2002a, Cornaglia & Balseiro 20026, Simon & 
Affleck 2002, Franzese et al. 2003, Cornaglia & Balseiro 2003, Kaul et al. 2005, Lewenkopf 
& Weidenmuller 2005, Kaul et al. 2006, Simon et al. 2006). 

One important recent development which has made very natural the idea that 
a magnetic impurity could be connected to finite size electron gas is the realization 
that Kondo physics was actually relevant to transport properties of quantum dots 
(Glazman & Pustilnik 2005). Indeed, as was pointed out by Glazman & Raikh (1988) 
and simultaneously Ng & Lee (1988), a quantum dot in the deep Coulomb blockade 
regime (so that particle number fluctuations are suppressed) containing an odd number 
of electrons, and sufficiently small that the temperature can be made much smaller 
than the mean level spacing between orbitals can be described by a Anderson impurity 
model, which, up to a Schrieffer- Wolff transformation, is essentially equivalent to a 
Kondo impurity. When the dot is weakly coupled to leads, these latter play the role of 
the electron gas. In the low temperature regime T <C a correlated state is formed 
which mixes the quantum dots and both drain and source wave-functions, leading to 
a large conductance (~ e 2 /h) in spite of the dot being in the deep Coulomb blockade 
regime. These predictions have been observed by Goldhaber-Gordon et al. (1998) a 
decade after these predictions, some later experimental realization even reaching the 
unitarity limit (van der Wiel et al. 2000). 

The great flexibility in the design and control of nanoscopic systems further 
motivated the study of many new configurations involving more exotic Kondo effects. 
One can cite for instance the theoretical design (Oreg & Goldhaber-Gordon 2003) 
and experimental observation (Potok et al. 2007) of the 2-channel Kondo, which 
has proved elusive in bulk system, the possible occurrence of SU(4) Kondo (Borda 
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et al. 2003, Le Hur & Simon 2003, Le Hur et al. 2004, Galpin et al. 2005, Hur 
et al. 2007) and its relevance to carbon nanotubes (Choi et al. 2005, Makarovski, 
Zhukov, Liu & Finkelstein 2007, Makarovski, Liu & Finkelstein 2007), or double dot 
system where Kondo physics might be in competition with RKKY interactions (Craig 
et al. 2004, Vavilov k Glazman 2005, Simon et al. 2005, Martins et al. 2006). 

The subject of Kondo physics and quantum dots is a very vast, and still rapidly 
developing, field. It is clearly not realistic to cover it in any reasonable way here, and 
I will only consider in more detail the, admittedly rather specific, aspects more closely 
related to the subject of this review. Indeed most of the more exotic designs imply 
at some point that a small quantum dot playing the role of a quantum impurity is 
connected to larger mesoscopic object, for which Kondo physics in itself is irrelevant, 
but such that finite size effects may become important. In other words, the context of 
Kondo physics and quantum dots makes it almost unavoidable to consider situations 
where the "quantum impurity" is connected to an electron gas for which finite size 
effects are important. 

As a consequence, for each such mesoscopic electron reservoir connected to the 
quantum impurity, two new energy scales enter into the description of the Kondo 
problem: the corresponding mean level spacing A R and Thouless energy E Th . The 
existence of a finite mean level spacing of the electron reservoir will, for instance, clearly 
modify the Kondo physics drastically for low temperature T <C Ar. This will affect 
the conductance (Thimm et al. 1999, Simon & Affleck 2002, Cornaglia & Balseiro 2003) 
as well as thermodynamic properties (Cornaglia & Balseiro 2002a, Franzese et al. 2003, 
Kaul et al. 2005), and considerable insight can be gained by considering the properties 
of the ground states and first few excited states of the system (Kaul et al. 2006, Kaul 
et al. 2008). 

The range of energy between A R and E Th is furthermore characterized by the 
presence of mesoscopic fluctuations (at all scales in this range) in the local density of 
states of the reservoir's electron, and thus by the fact that the density of states is not 
flat and featureless. In particular, one may wonder whether a (eventually fluctuating) 
Kondo temperature can be defined, and if physical quantities remain universal function 
of the ratio T/Tk- 

To fix the ideas, let us consider the s-d Hamiltonian (5.1) with a local density of 
states at the impurity site u(r = 0; e) = f £V |y?i(0)| 2 5(e — e«). In the semiclassical regime, 
u(r = 0; e) can, be written as the sum 

is(e) = u + i/ fl (e) (5.8) 

where uq is the bulk-like contribution (2.12) [one should of course include here either a 
realistic band dispersion relation or a cutoff at Do to account for the finite bandwidth], 
and the fluctuating term ffl(e) is a quantum correction associated with the interfering 
closed orbit contributions (2.15). Because ^fl(e) is the sum of rapidly oscillating terms, it 
will fluctuate not only with respect to the energy e, but also with respect to the position 
r of the quantum impurity or with respect to the variation of any external parameter 
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that may affect the classical actions Sj (2.4) on the scale h. As a consequence, one 
may think of ffl(e) as a statistical quantity with different realizations corresponding to 
various location of r or obtained by varying a external parameter in such a way that the 
classical dynamics remains unmodified, but that the phases exp(iSj/h) are randomised. 

Let us denote T^ the Kondo temperature of the associated bulk system for which 
u(e) is replaced by v> . For T it is possible to use a perturbative renormalization 

group approach in the same way as in the bulk, but including the mesoscopic fluctuations 
of the density of states. Following Zarand & Udvardi (1996), this gives in the one-loop 
approximation 

Jcff(Dcff) = i t r»»t / ^ ( ^ ' (5 ' 9) 

with 

■*<«> = ;£>;t^ ^ 

the temperature smoothed density of states (note the renormalization up to order two 
loops was given by Zarand & Udvardi (1996)). 

There are two different ways to use the above renormalization group equation. 
For some physical realizations, the fluctuations of the local density of states may be 
very significant, yielding even larger variation of the Kondo properties because of the 
exponential dependence in (5.4). In that case, one is mainly interested in the fluctuations 
of the Kondo temperature (which is now a functional of the local density of states) 
defined as the energy scale separating the weak and strong coupling regimes. One can 
then use the same approach as in the bulk and define T K \up\ as the temperature at 
which the one-loop effective interaction diverges, giving the implicit equation 



/•Do j, . 

J / —vp{u) = 1 . (5.11) 

Examples of systems for which the fluctuations of the density are large enough that 
one is mainly interested in the fluctuation of the scale T£ defined by (5.11) include, for 
instance, the case of "real" (chemical) impurities in a geometry such than one dimension 
is not much larger than the Fermi wavelength. This may be either a quantum point 
contact in a two dimensional electron gas (Zarand & Udvardi 1996) or a thin three 
dimensional film (Crepieux & Lacroix 2000). In that case there is a large probability 
that the impurity is located at a place at which the Friedel oscillations (i.e. the term 
in the oscillating part of the density of states associated to the trajectory bouncing 
back on the boundary and returning directly to its starting point) are large. In that 
case, the impurities located at a distance from the boundary such that the interference 
of the returning orbit are constructive have a significantly larger Kondo temperature. 
Disordered metals near or beyond the localisation transition provide another kind of 
systems with large fluctuations of the local density of states, and can lead in this 
particular case to a finite density of "free moments" for which the T£ [up] corresponding 
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to the space location is actually zero (Kettemann & Mucciolo 2007, Zhuravlev 
et al. 2007). 

In a typical ballistic or disordered-metallic mesoscopic system with all dimension 
much larger than the Fermi wavelength, the fluctuating part of the local density of states 
fa is a quantum correction to the secular term vq and is therefore parametrically smaller. 
As a consequence the fluctuations of T K [vp(u)] defined by (5.11) are not very large 
compared to T^-, which if it is taken only as an energy scale is somewhat meaningless. 

In the bulk case, Tk has however, beyond being an energy scale, another meaning 
which is to be the parameter entering into universal functions such as f x in (5.6). Within 
this framework, Tk is directly (and quantitatively) related to physical observables, and 
its fluctuations need not be large to be relevant. 

In the mesoscopic case, however, there is a priori no reason for any physical 
observable to be a universal function, since the number of "parameters" defining the 
problem can be considered as infinite (in the sense that one needs an infinite number 
of parameters to define the function v(e)). For a given temperature T > T\ however, 
it is possible to perform the perturbative renormalization group analysis leading to 
the effective interaction strength J cff (T), (5.9), which amounts to integrating out all 
degrees of freedom corresponding to energies larger than T. One can assume moreover 
that the features of the density of states at energies smaller than T will not affect 
the observables at temperature T. As a consequence, after renormalization, physical 
quantities can depend only on the parameters T, J e s(T), and i>p(T). Dimensional 
analysis then implies that dimensionless quantities such as Txi oc (T) can depend only 
on one single parameter, the product J{T) = vpiT) J eS (T). In particular Txi oc (^) 
should be equal for an arbitrary mesoscopic realization and for some bulk system with, 
at temperature T, the same value of the parameter J(T). This implies that if one 
defines the realization and temperature dependent Kondo temperature T K [vp]{T) by the 
generalisation of (5.7) to the mesoscopic case, namely 

* lm = ,mJt )M t)) ■ < 5 ' 12 > 

one obtains, within the one-loop approximation, a quantitative prediction for physical 
quantities in the mesoscopic case. For instance one has for the local susceptibility 

Txioc(T) = / X (T/TVM(T)) = / x (l/exp(J cff (TH(T)). (5.13) 

which relates the universal function f x , computed for the bulk to the impurity 
susceptibility in the mesoscopic case. 

From a practical point of view however, one may first note that with the 
definition (5.12), T^fz/^T^) = T£[i/g] is the same as the one obtained from (5.11). 
Furthermore, one finds that there is in general little difference between / x (T/Tx[^](T)) 
and f x (T '/T^fz/g]) in a large range of temperature above T^-. Therefore, although only 
the form (5.13) can be formally justified, I will continue the discussion assuming that 
in the high temperature regime T > T$-, the susceptibility can be described by the 
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universal form f x (T/T^) with a realization (but not temperature) dependent T£ defined 
by (5.11). 

It turns out that this approach gives a very precise prediction, provided Tj£ is 
computed from a fit or the two-loop approximation, and one uses T£ = + 57V with 
5TV defined by 

</o / — (v + 5vp(u;)) = J — v Q . (5.14) 

JT° K +ST K U JT° U 

Kaul et al. (2005) and Yoo et al. (2005) showed, by comparison with exact numerical 
quantum Monte Carlo calculations, that, even if one neglects the temperature variation 
of Tk[v/3}, the predictions from (5.13) are quantitatively extremely precise, even up 
to temperatures somewhat below T$-. This approach therefore provides a simple and 
quantitative way to discuss the mesoscopic fluctuations of the Kondo properties in the 
temperature regime T >T%. 

The nice thing here is that in this form, all the fluctuations of physical quantities 
in this regime are expressed in terms of the local density of states that we know how to 
relate to classical closed orbits. Indeed, (5.14) can be rewritten for 5T K <C as 

v / —v p (u) . (5.15) 



Furthermore using that vp{ui) = ui oc {eF + iuj) and the analytical continuation of (2.3) for 
complex energies in the same way as in (3.18), one can immediately relate the realization 
dependant Kondo temperature T^-[up\ to a sum over orbits starting and ending at the 
impurity site r = 



6Tk 1 ST n • 

t^ft — , — > Dj sin 



/ — exp(-tjuj/h) . (5.16) 



The last integral can be approximated by log ^-^J exp(— tjT^/H), thus providing a 

cutoff for trajectories with a time of return tj larger than h/T^. From this result, the 
variance of the Kondo temperature can then be obtained. Pairing a trajectory only with 
itself and, if the electron gas is time- reversal invariant, with its time-reversal symmetric, 
and using the sum rule (2.24), one has 

Tf ~ rrhvo J T% co[ co> J ^ l(t)e ' ^ 

with P^{t) the classical probability of return and /3rmt = 1 for a time-reversal non- 
invariant and 2 for a time-reversal invariant electron gas. 

For chaotic systems, using the model (2.26)-(2.27) for the classical probability of 
return then immediately gives 

1^1 K < 2 ^" ln2 >A . (5 . 18 ) 
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Interestingly, the same result can be derived from a random-matrix description 
(Kettemann 2004, Kaul et al. 2005). Indeed, using (5.15), the variance of 8T K can 
be expressed in terms of the correlator 

R^u^-J t ^4 M ^ M>) (5.19) 

= JL / / ^ v . (5 . 20) 

^iLL [</< - o 2 + m [(M-e«) 2 +^ 2 ] v 

Within the random-matrix model, the wave-functions have a Porter-Thomas 
distribution ((2.33) or (2.34)), and, neglecting the correlations between eigenstates, 
which introduce only 1/g corrections, one has ((\(p ai \ 2 — 1/ A) {\(p a2 \ 2 — 1/ A)) = 
8 ai a 2 (PRMT/ A 2 ) . Neglecting the fluctuations of the eigenergies e a and performing the 
integral, one then gets 

D / v /^RMtA f 1 iO X + L0 2 1 ^1+^2 

R 2 {uji,uj2) = — In 1 In 



7T \Ld 2 &i Wl oj 2 

and thus R 2 (T%,T%) = 2 ln2(/3 RMT A/7rT£), from which (5.18) derives immediately 
(Kaul et al. 2005). Within the random- matrix description, the full distribution of T K 
has been obtained by Kettemann (2004). 

To conclude this section, we see that the mesoscopic fluctuations of Kondo 
properties provide an example of physical systems where a relatively non trivial 
information is eventually encoded in the classical trajectories. Moreover through the 
sum rule (2.24) and the model (2.26)-(2.27) for the classical probability of return -P c e i(t), 
a connection with random-matrix theory can be made. 

It should be stressed however that the most interesting aspect of this mesoscopic 
Kondo problem, namely the fluctuations of physical properties in the deep Kondo regime 
(A <) T <C T^-, is still an open problem. This regime should be characterized by much 
larger fluctuations, which should be therefore easier to observe experimentally, as well, 
as has already been seen in exact numerical Monte Carlo calculations (Kaul et al. 2005), 
as a lack of universality. The perturbative renormalization group point of view used in 
the high temperature regime will clearly not be applicable at those temperatures, but a 
Fermi-liquid descriptions based on mean-field slave-boson techniques should nevertheless 
be able to shed some light on the physics dominating the mesoscopic fluctuations in this 
regime (Burdin 2007). 



6. Coulomb blockade peak spacing and ground-state spin of ballistic 
quantum dots 

I will turn now to a third (and last) illustration of physical systems for which the 
interplay between interference and interactions plays a dominant role, namely Coulomb 
blockade in ballistic quantum dots. Coulomb Blockade in itself is an essentially classical 
effect. It can for instance take place when a small metallic grain is weakly connected 
to electrical contacts and maintained at a temperature low compared to its charging 
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energy E c = e 2 /C (C is the total capacitance of the grain). In that case, in the lowest 
order in the grain-leads coupling, the conductance through the grain can be seen as 
a succession of transitions between states with different number of electrons N within 
the grain. Because of the Coulomb interaction between the electrons, these various 
states will have an energy difference of order E c , and in general the conservation of the 
total energy cannot be fulfilled. As a consequence, the transport through the grain is 
blocked. If, however, the grain is capacitively coupled to an external gate, this latter 
will also affect in a different way the energy of the states with a different number of 
electrons. The potential of the gate, V g , can therefore be tuned so as to adjust the 
energies of the various states. Thus, as a function of V g , the conductance through 
the grain will display an alternation of valleys and peaks (for reviews see (Grabert & 
H. Devoret 1992, Kouwenhoven et al. 1997)). 

The Coulomb-blockade process described above does not imply a quantum 
mechanical effect, and can in particular be observed if the mean- level spacing, or even 
the Thouless energy, is much smaller than the temperature. Using very small systems, 
such as ballistic quantum dots built in GaAs/AlGaAs, one can nevertheless increase 
sufficiently the one-particle mean-level spacing A so that A > T, and the transport in the 
dot takes place through a single (or a few) levels. In this regime, one observe fluctuations 
in both spacings and heights of the conductance peaks, fluctuations which encode non- 
trivial information about the many-body ground-states (with various particle numbers), 
and possibly a few excited states. 

6.1. Constant-interaction model and experimental distributions 

As far as the peak-height fluctuations are concerned, predictions (Jalabert et al. 1992) 
based on the so called "constant-interaction model" appeared to fit the distributions 
measured in the earlier set of experiments (Chang et al. 1996, Folk et al. 1996). 
In this model, beyond a classical charging term (Ne) 2 /2C , interactions among the 
electrons are completely neglected, and Porter- Thomas (i.e. RMT-like) fluctuations are 
assumed for the one-particle wave functions. Even at this early stage, however, the 
presence of correlations between successive peak heights appeared incompatible with a 
simple RMT description, and pointed to the role of short periodic orbits (Narimanov 
et al. 1999, Kaplan 2000, Narimanov et al. 2001). Furthermore, more recent experiments 
(Patel et al. 1998) showed decreased probability of having either very large or very small 
heights. Several suggestions have been put forward to interpret these deviations, among 
which the effect of inelastic scattering (Rupp et al. 2002), of spin orbit (Held et al. 2003), 
or of interactions (Usaj & Baranger 2003). It remains that for peak-height distributions, 
the constant-interaction model seems to capture a good part of the relevant physics. 

The status of the peak-spacing fluctuations, however, is very different. For this 
quantity, the constant-interaction model gives very striking predictions, since, because 
of the spin degeneracy, the alternation between singly and doubly occupied orbitals, as 
the number of electrons N in the dot increases, is associated with a strongly bimodal 
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distribution (Sivan et al. 1996). Indeed, in the zero temperature limit, the position V* 
of a conduction peak is determined by the energy conservation condition 

S N (V g *)+v = £ N+1 (V g *) (6.1) 

(n is the chemical potential in the leads). Writing S N (V g ), the ground-state energy of the 
dot with N electrons, as the sum of some "intrinsic" part plus a term — (C g /C)eNV g 
due to the capacitive coupling with the control gate (C g is the capacitance of the gate to 
the dot), one obtains that the spacing in V g between two successive peaks is proportional 
to the second (discrete) energy difference 

(V g *) N ^ N+1 - (V g *) N _^ N oc 5 2 E^ ^ E^ +1 + Et 1 - 2E» . (6.2) 

In the constant-interaction model, the ground-state energy of the dot with N 
electrons is written as 

occupied icr 

where the last term is simply the one-particle energy of a system of non-interacting 
fermions (i and a are respectively the orbital and spin index and q the corresponding 
one-particle energy). One obvious prediction of this model is that the ground state spin 
of the dot can be only (for N even) or 1/2 (for N odd). Furthermore, because of the 
spin degeneracy, one gets 

5 2 E* = e 2 /C for odd N , (6.4) 

5 2 E* = e 2 /C + (e N/2+ i - ejv/ 2 )for even N , (6.5) 

and the peak spacing is the superposition of a Dirac delta function (corresponding to odd 
N, referred below as "odd spacings") and of the nearest neighbor distribution P mis (s) 
of the 6i (corresponding to even N, referred below as "even spacings"). For a chaotic 
system, this is described by random-matrix theory, giving (2.32). 

The distributions observed experimentally share very little resemblance with 
this prediction. The first set of experimental results (Sivan et al. 1996, Simmel 
et al. 1997, Simmel et al. 1999) not only did not show any trace of bimodality, but 
the width of the distributions seemed on a scale of the charging energy rather than on 
the one of the mean- level spacing (as P nns (s) was expected to be). It was soon realized 
that these very large distributions were dominated by switching events, that is, by the 
displacement of trapped charges located between the quantum dot and the control gate. 
Further experiments (Patel et al. 1998, Liischer et al. 2001, Ong et al. 2001) where then 
performed which took care of maintaining the level of the noise associated to switching 
events well below the mean-level spacing. 

In the work of Patel et al. (1998), seven dots of different sizes (and two different 
densities) have been investigated for finite as well as zero magnetic field. Measures were 
taken at an electron temperature of 100 mK, which corresponded to a ratio T/A ranging 
from 20% for the smallest dot (largest A) to 60% for the largest dot (smallest A). The 
noise level, estimated by comparing the data when reversing the sign of the magnetic 
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Figure 7. Histogram of the peak-spacing data for the dot 1 of Patel et al. (1998). The 
noise level and temperature correspond respectively to 8% and 20% of the mean-level 
spacing A. The electronic density n s = 2. 10 11 cm -2 correspond to a gas parameter 
r s ~ 1.24. This figure is taken from Patel et al. (1998). The dotted and solid lines 
corresponds to predictions from the constant interaction plus spin resolved RMT model 
(CI+SRRMT) presented in this reference. 

field, ranged from 8% to 94% of the mean-level spacing. The peak-spacing distribution 
obtained for the smallest dot, which turns out to be the one showing the less noise 
(T/A = 20%, noise = 8% of A), is reproduced in figure 7. As seen there, this distribution 
remains essentially uni-modal and Gaussian, and therefore clearly incompatible with the 
constant-interaction model predictions. 

Data from the same series of dots have been re-analysed by Ong et al. (2001). In 
this latter work, special emphasis has been put on the study of the difference between 
odd and even spacings predicted by the constant-interaction model. Some odd-even 
effect could actually been observed for the dots exhibiting the least amount of noise 
and the smallest T/A ratio. This is illustrated in figure 8 which shows separately the 
distributions of odd and even peak-spacings for the same dot as in figure 7, but with 
a slightly different gate configuration so that the effective area and thus the ratio T/A 
are somewhat larger (T ~ 30% A). In that case, some differences are seen between 
the odd and even distributions. They remain however clearly incompatible with the 
constant-interaction model predictions (6.4)-(6.5) . In particular, the total peak-spacing 
distribution that would be obtained by summing the two curves would, as the one shown 
in figure 7, be uni-modal and essentially Gaussian. 

The same kind of distributions were obtained by Liischer et al. (2001) for a dot 
defined by local oxidation with an atomic-force microscope (rather than electrostatic 
gating as in Patel et al. (1998)). The steepest confining potential obtained in this way 
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Figure 8. Probability distributions of even (solid line) and odd (dash line) peak 
spacings. The data correspond to the dot 1 of Patel et al. (1998), but for a different gate 
configuration which is such that, in spite of the slightly lower temperature T — 90mfc, 
the ratio T/A is about 30%. Some visible difference is seen between the odd and even 
distributions, but these latter clearly do not correspond to the constant-interaction 
model predictions (6.4)-(6.5). This figure is taken from Ong et al. (2001). 

made it possible: i) to limit the shape variation of the dot as the lateral gate voltage V g 
is changed; ii) to use a back gate to increase the electronic density (n s ~ 5.9 10 11 cm~ 2 ) 
and thus reduce the gas parameter to a value r s ~ 0.72; and furthermore to define a 
smaller quantum dot so that the electronic temperature T = 120 mK used amounts 
to a significantly smaller ratio T/A ~ 5%. As illustrated in figure 9 the peak-spacing 
distributions nevertheless remained uni-modal and essentially Gaussian. 

Other experimental evidences that the constant-interaction model does not properly 
account for the ground-state properties of quantum dots can be obtained by considering 
their total spin. Indeed, as already mentioned, this model implies a "naive" occupation 
of the orbitals, and therefore, for the ground state, a total spin zero for even or one half 
for odd number of particles respectively. Already in Luscher et al. (2001), the parametric 
variation of the Coulomb-blockade peak position as a function of a weak perpendicular 
magnetic field, and in particular the existence of kinks (Baranger et al. 2000), pointed 
to the existence of ground-states with spin one. Such spin-one ground-state were further 
used to observe Kondo effect in small dots with even number of particles (Kogan 
et al. 2003). More systematic studies of ground-state spins were later performed using 
stronger in-plane magnetic field coupled to the spin of the electrons through Zeeman 
effect (Rokhinson et al. 2001, Folk et al. 2001, Lindemann et al. 2002). Not enough data 
have been collected in this way to obtain an experimental distribution of ground-state 
spins. Folk et al. (2001) have nevertheless obtained clear evidence of the existence of 
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Figure 9. Histogram of the peak-spacing distribution for B = (a) and B ^ (b) 
for the quantum dot studied by Liischer et al. (2001). The various solid curves are fit 
to theoretical predictions discussed in this reference. Inset: the histogram is the same 
as (b), but with a different scale. This figure is taken from Liischer et al. (2001). 

non-trivial (i.e. different from zero or one half) ground-state spins for semiconductor 
quantum dots similar to the ones used for figures 7 and 8. 

6.2. The universal Hamiltonian 

These discrepancies between the constant-interaction model predictions and what was 
observed pointed clearly to the role played by the residual interactions (i.e. beyond 
simple charging energy) in these systems, and motivated a large number of works ranging 
from exact numerical calculations on small system (Sivan et al. 1996, Berkovits 1998) 
to various kind of self-consistent approximations :(:. In the end however, it seems that, 
even if not all aspects of the experimental data can yet be explained, the main features 
are compatible with a relatively simple Fermi-liquid description provided the role of the 
spins are properly taken into account (Blanter et al. 1997, Brouwer et al. 1999, Baranger 
et al. 2000, Ullmo & Baranger 2001, Aleiner et al. 2002). Once the consequences of a 
finite temperature are included (Usaj & Baranger 2001, Usaj & Baranger 2002), the 
predictions for the distributions of Coulomb-blockade peak-spacing are actually in very 
decent agreement with at least some of the experimental data. 

What is meant by "Fermi-liquid" description here is the picture given in section 2.3: 

| As for instance Hartree Fock (Walker, Gefen & Montambaux 1999, Walker, Montambaux & Gefen 
1999, Cohen et al. 1999, Ahn et al. 1999) or density functional theory (Stopa 1993, Stopa 1996, Koskinen 
et al. 1997, Lee et al. 1998, Hirose & Wingreen 1999, Hirose & Wingreen 2002, Jiang et al. 20036, Jiang 
et al. 2003 a). 
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an effective "mean field" one-particle Hamiltonian H^f which, for a specific system 
could be in principle obtained by minimising the Thomas-Fermi functional (2.48) (see 
for instance the discussions given in Ullmo et al. (2001) and Ullmo et al. (2004)), and a 
weak residual interaction well approximated by the RPA-screened Coulomb interaction, 
or even its long wave length zero frequency limit (2.46)-(2.47). As pointed out in 
section 2.3 there is up to now, for a generic system, no real derivation of this picture 
from basic principles. The point here is thus not to define an approximation scheme 
to predict quantitatively the energies of a few specific realizations, but rather to have 
a model which, once supplemented by some modelling for the fluctuations of the one- 
particle eigenstates <fi and energies typically RMT for chaotic systems, gives correct 
predictions for statistical properties such as the peak spacing distribution, etc.. 

6.2.1. Time-reversal non-invariant systems As for orbital magnetism (see section 3.2), 
time reversal invariance introduces a slight complication in the discussion because of the 
Cooper series for which, at r s ~ 1, high-order terms need to be included and eventually 
renormalise the effects of the interactions. I shall therefore start the discussion assuming 
the presence of a magnetic field strong enough to break time reversal invariance, but 
nevertheless small enough not to change qualitatively the classical dynamics within the 
system (and in particular so that it remains in the chaotic regime). In that case, the 
screened interaction can be treated at the first order of the perturbation. 

In the limit where the residual interactions are neglected, the many-body eigenstates 
are Slater determinants characterized by the occupation numbers n ia = 0, 1 of the mean- 
field Hamiltonian's wave-functions <fi . For the states such that all singly occupied 
orbital have the same spin polarisation (which as we shall see are the only ones that 
can be the ground-state of the system), the usual non-degenerate perturbation theory 
can be used: in first order the eigenstates are unmodified and have an energy (Ullmo & 
Baranger 2001) 

E{n ic ) = E sm (N) + rua* + E m {n m } (6.6) 

icr 

E K1 {n ia } = ^ nurMijUja* - ^ ^ n ia N i:j n ja . (6.7) 

ia,ja' i,j,a 

Here E sm (N) is a smooth contribution containing essentially the electrostatic energy 
(Ne) 2 /2C , so that the first two terms of (6.6) correspond to the constant-interaction 
model. The residual interaction term E m {n ia } is then expressed in terms of 

M tJ = J drdr'\^{r)\ 2 V sc {r - r')\ Vj {r')\ 2 , (6.8a) 

N tj = J drdr' Vl (r) V *(r)V sc (r - r>,(r>*(r') . (6.86) 

We shall see below (see section 6.3) that the fluctuations of the M^- and Ny are 
parametrically smaller than their mean values (M^) and (N^-). We shall therefore 
for a moment neglect the former. The quantities (M^) and (N^) can be derived 
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by different methods, in particular from the knowledge of the two-point correlation 
function (2.38a) and a Gaussian hypothesis for the higher-order correlation functions 
(Srednicki 1996, Hortikar & Srednicki 1998). This can also be done directly from the 
random-plane-wave model. For instance, for a billiard system with only a kinetic energy 
term in a region of space of volume (or area) A, the insertion of (2.39) in (6.8a) gives 

{Mij) = J drdr' ( a in a in> a jv a * v >) ex P Q[ r (*V ~ *V) + r '(Pj> ~ P^')Vj ■ ( 6 - 9 ) 

The Gaussian character of the a iit and (2.40) imply (ai^a* a^a^) = (a^a* ,) (a^a*^,) + 
Sij(ai^a* iri ,)(a iri a*^) = (AM)' 1 (8 ^8^ + 8^8 ^8^). One obtains in this way 

(Mij) = (M#j) + 6ij 8(M U ) , 

(M^> = iv sc (0) (6.10) 
<f<M«) = -^(V; c )f. s . (6.11) 

where (t^ c )f. s . is the average of V sc (q) for q = p—p', p and p' being uniformly distributed 
on the Fermi surface. For two and three dimensions: 

(V; c )f.s. = t^ J d9V sc (A; FV /2(l + cos#)) (d=2) , (6.12) 

= ^- J sm(9)d9dipV sc (k F ^2(l + cos9^j (d=3) . (6.13) 

A direct calculation gives similarly 

(N i ^-> = <J(M ii ) = i(y BC ) f . s .. (6.14) 

The equality between (Nj^j) and 8{M U ), which may seem somewhat surprising at first 
sight, can be understood as taking root in the invariance of the residual-interaction 
energy under a change in the spin quantization axis. Finally Njj = Ma and thus 
compensates the corresponding a = a' contribution. 

Using the approximation (2.46)-(2.47) for the screened interaction one gets 

^sc(O) (6.15) 

gs 

^(V sc )i,. = JrpaA, (6.16) 

where A is the mean- level spacing between the q's and g s = 2 the spin degeneracy. 
The last equation actually defines the dimensionless constant Jrpa, which, according 
to (2.46)-(2.47), is seen to be of order r s in the high density (r s — > 0) limit but of 
order one (although always smaller than 0.5) as r s becomes of order one. More properly 
however, (t^ c )f. s ./A should be interpreted as the parameter /° of the Fermi-liquid theory 
(Pines & Nozieres 1966). This ratio has been computed, as a function of r s , in a variety 
of ways. Figure 10 shows for instance, for two dimensional systems, a comparison 
between the value Jrpa obtained from the RPA (actually Thomas-Fermi) expression 
(2.46) and the numerical Monte-Carlo evaluation from Tanatar & Ceperley (1989). For 
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Figure 10. Comparison between the value Jrpa obtained from the RPA (actually 
Thomas-Fermi) expression (2.46) and the numerical Monte-Carlo evaluation of /° from 
Tanatar & Ceperley (1989). (Courtesy of Gonzalo Usaj.) 



r s — 1 for instance, we see that the Monte-Carlo result ~ 0.34 is slightly larger than 
(V sc )f. a ./A = 0.28 obtained from RPA. 

A few remarks are in order here. First, we note that, as expected, the mean value of 
the interaction contribution to the total energy is much smaller than the mean value of 
the non-interacting contributions (whether the electrostatic term E sm (N) ~ (eN) 2 /2C 
or the smooth part of the one-particle energy ^2n ia ei). In the same way, we shall see 
below that the fluctuations of the interaction contributions are parametrically smaller 
than that of the non-interacting term. However, the mean values of the My's and N^-'s 
are on the same scale (A) as the fluctuations of the one-particle energies 6j. In the end, 
it is the interplay between these fluctuations of a large term, and the mean value of a 
smaller one which makes possible a qualitative change from the interactions although 
they can be treated perturbatively. 

The second point that should be stressed is that, since the variation with the 
orbital indices of the Ny's and M^-'s is neglected, the residual-interaction term (6.7) 
depends only on the number of same-spin and different-spin pairs of particles, which 
themselves can be expressed in terms of the total number of particle N and of the 
difference (N + —N_) between the numbers of majority and minority spins. We consider 
here only states such that all singly occupied orbitals have the same spin polarisation, 
hence (iV + — AL) is just twice the total spin S of the system. Keeping in mind that 
5 (Ma) = (Njj), we get from simple algebra 



E m (N, S) 



N(N-l) 



N(N — 4) 



4 



(%)/2-5(5+l)(N ij ) . 



(6.17) 
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The A-dependent part can be safely aggregated with the large smooth term E sm (N), 
and it turns our that the residual-interaction term is in this case proportional to the 
eigenvalue S(S + 1) of the total spin square S 2 ot of the dot. 

As was pointed out by Aleiner and coworkers (Kurland et al. 2000, Aleiner 
et al. 2002) this result is to be expected, and is just the consequence of the general 
symmetries of our problem. Indeed, the description we have used for the eigenstates 
of the non-interacting system is, within an energy band centered at the Fermi energy 
and of width the Thouless energy, equivalent to a random-matrix model. This means 
that as long as one considers an energy scale smaller than the Thouless energy (which 
is assumed here), there is no preferred direction in the Hilbert space, and any physical 
quantity should, after averaging over the ensemble, yield a result which is invariant 
under rotation of the Hilbert space. As a consequence, if fluctuations are neglected, the 
residual-interaction Hamiltonian 

^m = 9^ VijMclc^c^Ckv (6.18) 



i,j,k,l 
aa' 



with 



V^ki = I drdr'(p*(r)(p*(r')(pi(r')ip k (r) 



(6.19) 



should have a "universal form" , i.e. should be expressed only in terms of the invariants, 
which, without time-reversal symmetry, are the number operator A and the total spin 
square S 2 



tot- 



Using the random-plane-wave model we can get an explicit expression of this 
universal-Hamiltonian form, including the values of the various parameters. Indeed, 
inserting (2.39) in (6.19) one has 

Vi jk i = a *i^ a *ju a iy' a k^ j drdr'exp ( ^[r • (jy ~ Pv) + r ' ' (fV ~ Pv)\ J V sc (r-r') .(6.20) 



With (a*^a* ju a lv ,a kfM ,) = (a*^a k ^)(a* ju a lu ,) + (a*^a lu ,)(a* ju a kfM ,) = ^[8 ik 8^ -<^<W+<^<V 
8jk8 U u'}, one gets 



aa' 



Vsc(0) 

A 



C i<T C ja ,C j°~' C i°~ 



(V sc ) 



Ls -J J 



C icr C j<7 ,C icr' C jcr 



We can then use the equalities 

EiSv c ^'Q CT = A>(A>-l) 



1,3 
aa' 



E c l c l' c ^ = - + 2N - 2S 



2 

tot 



1,3 
aa' 



+ fluctuating terms .(6.21) 



to write the residual-interaction part of the Hamiltonian as the sum of two terms 



Hri = H 



(AT) 
RI 



+ H$. The A-dependent part = (U SC (0)/2^)A(A - 1) 



(TV) 
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((V SC ){. S ./A)[(N /4)(iV - 4)] can be aggregated with the smooth term E sm (N) in (6.6), 
while the spin part can be expressed as 

Hg> = -^^S t 2 ot = -A /o «S t 2 ot . (6.22) 

We obtain in this relatively pedestrian way that, as expected, when fluctuations are 
neglected, from general symmetry considerations, the residual interaction is, up to a 
term which depend smoothly on the total number of particles, proportional to Sf ot . 
We moreover obtain that the proportionality constant is — (V sc )f. s JA which should, as 
before, be interpreted as —A fg. 

We now understand that the basic ingredient missing in the constant-interaction 

* (s) 

model was the existence of the term H^j . This term tends to polarise the spin of 
the quantum dot, and is of the same order (the mean-level spacing A) as the cost in 
one-particle energy associated with moving a particle in higher orbitals. Because the 
numerical value of /° is smaller than one half, the one-particle energy cost is on average, 
always numerically larger than the spin polarisation term. However the fluctuations of 
the one-particle energies taking place also on the scale A, imply that, depending on 
the actual value of the q's, the ground state may actually have a non-minimal spin. This 
in return will modify the quantities such as the peak-spacing distributions. Figure 11 
shows the distributions obtained within this model for r s — 1. 



6.2.2. Time-reversal invariant systems Physical effects not contained in the universal 
Hamiltonian are going to further modify the peak-spacing and spin distributions. 
Before considering them, let us return to the question of time reversal invariance, 
assuming again chaotic dynamics. When this symmetry is present the random matrix 
is no more invariant under all unitary transformations, but instead under the smaller 
group of orthogonal transformations. In addition to N and S, such transformations 
leave also T = X]j c «T c «l invariant. Taking into account that the residual-interaction 
Hamiltonian should be invariant under a global phase change of the unperturbed one- 
particle eigenstates (and thus of the implying that only the product T^T should be 
involved), and that H mt contains only four-fermion products, we see that an extra term 

iZg? = -J T f ] f (6.23) 

should be added to the general form of the universal Hamiltonian (Kurland et al. 2000, 
Aleiner et al. 2002). 

Using the form of the plane-wave model valid for time-reversal invariant systems 
[see discussion below (2.40)], we can as before express the interaction matrix element 
V ijk i through (6.20). Now, however, (a^a^a^a^) = ^[5^5^ ■ 5 jt 5 uv/ + SuS^ ■ 
^jk^un'+^ij^n-v-^ki^u'-n'], where the last contribution is specifically due to time-reversal 
invariance. One obtains 

(Vnjj/TRI = (Viiii)TRmO~ii + —^{V 8C )t mSm , (6.24) 

where (V^^trni is the mean value of the completely diagonal matrix element in the 
time-reversal non-invariant case. 
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Figure 11. The probability density of Coulomb-blockade peak spacings obtained from 
the universal-Hamiltonian model at zero temperature. The total distribution (solid) 
as well as that for N even (dashed) and odd (dotted) is given for two strengths of 
the interactions: (a) r s = 1 and (b) the zero-range interaction limit V sc (r — r') = 
g s _1 ^o — r ')- The presence of a 5-function in the distribution is particularly clear 
in the cumulative distribution functions (the integral of the probability density) shown 
in (c). Insets show the probability of occurrence of ground-state spins in the two cases. 
The spacing is in units of the mean- level separation A, and the origin corresponds to 
the classical spacing e 2 /C. This figure is taken from Ullmo & Baranger (2001). 
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We obtain in this way that Jt equals J5. What makes somewhat more difficult the 
discussion of systems with time-reversal invariance is the need to consider higher-order 
perturbation terms. Indeed, while N and S 2 were commuting with the unperturbed 
Hamiltonian, the product T^T is not. In particular, the evaluation of second-order 
corrections in the residual interactions (see (6.27) in the next subsection) shows that 
the matrix elements associated with the promotion to an empty orbital of two electrons 
(with opposite spins) occupying the same orbital are now of the order of the mean-level 
spacing A. This correction is thus by no means small. 

We recognise however in H^p the usual pairing Hamiltonian used in the study 
of superconductivity (except that the interaction here is repulsive). It is well known 
(Abrikosov et al. 1963) that the important higher-order terms are, as in section 3.2, the 
Cooper series shown in figure 1. In the same way as for the magnetic response, the main 
role of these higher-order terms is therefore to renormalise the interaction in the Cooper 
channel according to (4.44) (for d = 2). This argument is often used to neglect 
since the renormalised coupling constant goes to zero in the limit g — > 00. However, 
for typical numbers of particles in the studied experimental dots (e.g. from 340 to 1000 
for Patel et al. (1998) or about 200 for Liischer et al. (2001)) the "large" logarithm 
g s _1 In (k F L) ~ g s _1 In (47riV/g s ) remains in the range [3.5, 4.4]. The higher-order terms 
therefore reduce the effects of the interactions in the Cooper channel (i.e. a smaller 
effective parameter Jt should be used) but do not eliminate them completely. 

The role of the renormalization in the Cooper channel is presumably the reason 
why density functional calculations in the local spin density approximation, while they 
agree very well with the kind of Fermi-liquid description given here for time-reversal 
non- invariant systems (Ullmo et al. 2005), significantly overestimate the role of the 
residual interactions in the time-reversal invariant case. The higher-order terms appear 
not to be treated correctly in this approach (Jiang et al. 2003 b, Jiang et al. 2003a, Ullmo 
et al. 2004). 

6. 3. Beyond the universal Hamiltonian 

In figure 11 one sees that predictions derived at zero temperature from the universal 
Hamiltonian are noticeably different from the ones ((6.4)-(6.5)) obtained within the 
constant-interaction model. In particular, as shown in the inset, the exchange term 
(6.22) may give rise to a non-zero proportion of non-naive (i.e. different from zero or 
one half) ground-state spins. However, even if one accounts for some experimental noise 
that would smooth some of the sharpest features, the peak-spacing distributions derived 
in this way do not have the uni-modal / Gaussian-like shape observed experimentally. 
It is clear that the universal Hamiltonian cannot be the full story. Nevertheless it 
seems not unreasonable to assume that it provides a decent starting point, which needs 
to be supplemented by a few other physical effects before a complete description is 
reached. Among the possible candidates, come naturally to mind: i) scrambling (Blanter 
et al. 1997, Ullmo & Baranger 2001, Usaj & Baranger 2002, Jiang et al. 2005) and gate 
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effects (Vallejos et al. 1998, Ullmo & Baranger 2001, Usaj & Baranger 2002, Jiang 
et al. 2005); ii) fluctuations of the residual-interaction terms (Blanter et al. 1997, Ullmo 
& Baranger 2001, Usaj & Baranger 2002), and contributions beyond order one (Jacquod 
& Douglas Stone 2000, Jacquod & Stone 2001, Usaj & Baranger 2002) ; and finally iii) 
temperature effects (Usaj & Baranger 2001, Usaj & Baranger 2002). I shall consider 
each of them separately, and end with a discussion about the importance of the more 
or less chaotic character of the actual dynamics within the dots (Ullmo et al. 2003). 

6.3.1. Scrambling and gate effects In addition to residual- interaction terms (6.7), there 
is one noticeable difference between the Fermi-liquid approach that we follow here and 
the constant-interaction model. Indeed, we have in principle a well defined procedure 
to specify the self-consistent confining potential U m f(r) seen by the electrons: for a 
given experimental configuration, it can be obtained by solving the Thomas-Fermi 
equations (2.51)-(2.52). [Note it is important not to use here for U m f(r) the potential 
obtained from a self-consistent Hartree Fock, or density functional, calculation, as this 
would mix the scrambling and the gate effects with (some treatment of) the residual 
interactions]. As one extra electron is added into the quantum dot, or as the voltage of 
the control gate is changed (without changing the number of electrons in the dot) from 
V*(N - 1 -> N) to V*(N — > iV + 1), the potential U m{ (r) will be modified to a new 
value U m f(r)+8U m f(r). As a consequence, the one-particle energies ej, and therefore the 
corresponding contribution to the ground-state energy, will be modified. The fluctuating 
part of the this energy change is referred to as scrambling (for the variation of N) and 
gate effect (Vallejos et al. 1998) (for the variation of V g ). 

An evaluation of the magnitude is possible analytically for the simple geometry of 
a circular billiard, with in addition (for the gate effect) the assumption that one is using 
an "universal gate" , i.e. one which is featureless and of the same size (or larger) than the 
quantum dot. Under these assumptions, both effects appear to have a variance scaling 
as 1/g, to be of the same size, and to be comparable in magnitude with the fluctuations 
of the residual interactions that I shall discuss below (Ullmo & Baranger 2001, Usaj & 
Baranger 2002, Aleiner et al. 2002). 

A more careful analysis reveals that a realistic modelling of the geometry of the 
quantum dots, which implies the numerical solution of the Thomas- Fermi self-consistent 
equations, gives a drastically different picture (Jiang et al. 2005). Indeed because in 
practice lateral plunger gates are used rather than the kind of "universal" ones assumed 
above, the quantum dots are affected asymmetrically by the gates, and gate effects turns 
out to be stronger than scrambling. On the other hand, the very smooth confinement 
associated with a realistic potential U m f(r) is much easier screened than the hard- wall 
boundaries of a billiard. This turns out to diminish noticeably the magnitude of both 
scrambling and gate effects: In the configuration studied in Jiang et al. (2005), very 
small fluctuations (less than a percent of the mean-level spacing for the root mean 
square (r.m.s.)) is obtained for the scrambling, so that it can presumably be neglected 
in actual experiments. Because of the use of a lateral plunger gate, gate effects are - 
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again in the configuration studied in Jiang et al. (2005) - not as small as scrambling, 
giving rise to fluctuations of the order of 7% of a mean-level spacing. This figure remains 
however smaller than the noise level for all the dots studied in Patel et al. (1998) (though 
marginally smaller for the quieter one, used for figure 7). 

Since some particular care has been exercised by Liischer et al. to minimise gate 
effects (the lateral gate is not so narrow, and the higher electronic density makes 
screening more effective) it can be assumed that the quantum dot studied in Liischer 
et al. (2001) shows even smaller gate effects than their evaluation in Jiang et al. (2005). 



h vPm + Pu ~ Pt ~ Pr,) {r - r') 



.(6. 



6.3.2. Fluctuations of the residual-interaction terms and higher-order correction 
Screening and gate effects correspond to changes of the one-particle energies when the 
smooth effective confining potential is modified for one or another reason. They lead 
to variations of the parameters of the universal Hamiltonian, but remain in some sense 
contained within this description. 

One may, however, consider terms, associated with the fluctuations of the residual 
interactions, which are corrections to the universal Hamiltonian. For instance, the off- 
diagonal (i.e. 7^ (k,l)) coefficients Vijki (6.19) are zero on average, but have a non- 
zero variance which can be evaluated easily within the random-plane-wave model. For 
instance if % ^ k, I and j ^ k, I, considering for simplicity the zero range approximation 
V sc (r — r') = t / (7 1 5(T — r'), and using (2.39) and (2.40) together with the Gaussian 
character of the a^'s, one can write 

11 f ' 

(i ^ |2)= ^(M^^y^ r/exp ' 

Once the momenta are constrained to be on the Fermi circle (or sphere, this seems to be 
irrelevant here), the condition p^ + p v — p T + p T] implied by the integration over space 
imposes p^ = p T + 0(h/L) and p u = p v + 0(h/L) or the converse (the 0(h/L) originates 
from the width given to the Fermi surface in the RPW model). As a consequence, using 
that the number of plane waves M scales as the dimensionless conductance g, we get 

oc 1 . (6.26) 
A z g z 

We see that the smallness of the typical size of the off-diagonal elements is not so much 
due to the weakness of the interactions (in some sense, UQ X 5(r — r') is of order one once 
put in the proper units), but rather originates from the self averaging associated with 
the integration over space of fluctuating quantities. 

As the matrix elements are small in the semiclassical limit g — > oo, one can expect 
for instance that second-order correction to the ground-state energy 

E m v l(*fl^K)l 2 (6 27) 

E ° ~ 1^ F (o) F ( 0) ( 6 - 27 ) 

also scales as some power of 1/g. The analysis is made somewhat more complicated by 
the fact that the number of one-particle levels within the Thouless energy grows with g, 
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and thus also the number of one-particle-one-hole and two-particle-two-hole excitations 
involved in the summation. A careful analysis performed by Usaj & Baranger (2002) 
shows nevertheless that the typical size of the second-order correction to the addition 
energy scales as A/g, is numerically already extremely small (less than 2% of A) for 
N = 500, and can be safely neglected. This furthermore gives us extra confidence that 
the kind of Fermi-liquid perturbative approach that we are following is indeed valid. 

The largest terms neglected by the universal-Hamiltonian description is therefore 
associated with the fluctuations of the diagonal residual interaction terms = 
and = V^j that we are going to evaluate more carefully now in the case of a two 
dimensional quantum dot. 

As before, we shall base our calculation on the random-plane-wave model. It turns 
out that the variances of the M^-'s and N^-'s (expressed in units of the mean- level 
spacing A) keep a dependence on the size of the system. To reproduce correctly this 
size dependence one needs the second version of the random-plane-wave model, in which 
the condition (2.41) is imposed on the states used in the expansion (2.39), giving a width 
~ h/L to the Fermi surface. In this way, one can write the My's as 

Mjj = a «m%J2%M3S>4^PMl-PM2>-PM 3 +PM4^ sc (PMl ~~ ) 

Ml,^2,/i3,£t4 

= X^MW^W? (6.28) 

Q 

where one introduces the definition 

W ig = £ tvx^- (6-29) 

Pmi-Pm 2 =< J 

The random-matrix model implies 

(aifeiO* fe2 ) = {k F L)~ 1 8 i j8p l p 2 if 5k < tt/L (6.30) 
= if 5k > tt/L 



with 5k — | ] fci| — ki\ and ki = y/2mei/H. From this, one deduces 

var(M^) ^^J2^ 2 M(m q \ 2 )(\W 3q \ 2 } . (6.31) 

|Wjq| 2 can be interpreted as (2irki)~ 2 times the area of the intersection of two rings of 
diameter ki and width 2n/L, centered at a distance \q\ = q. Simple geometry gives, for 
2tt/L < \q\ < 2ki - 2tt/L 

mq]2) ~ ( g L)V(2W-(3p ' (6 ' 32) 
One obtains for % ~ j 

™ {U '*^^U i mv=m - (6 ' 33) 
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The variance of Mjj and and the covariance between M^- and can be 



computed along the same lines, and one gets 

2 f 2k ^' L dq [V sc (q) + V sc (2k)] 2 



var(N^) ^ -5L ^ 



/L q (2kLy-(qL) 



2 



(6.34) 



/iwr ivr \ /m wn w 4 /" 2fc - VL ^ c (g)fe(g) + v sc (2^ 

(M^-N^) - <M^)(N^> _ _ - {2kL)2 _ {qL)2 . (6.35) 

The diagonal part of the direct residual interaction has an extra contribution because 
of the additional fashion in which the wave-functions can be paired: 

va. r( M„) * M^) + -p I f (2^)^ (^ ' < &36 » 

In the zero-range interaction limit, the expressions for the variance of the M's and TV's 
simplify considerably and one finds in this case 

var(M y ) = var(N ¥i ) = ^^^(1 + • (6-37) 

Note that the decay of the wave-function correlations appearing in the variance produces 
a factor of l/k F L in the root mean square compared to the mean. The In(kL) factor is 
special to two dimensions; it comes from the 1/kL decay of the wave-function correlator 
in this case. 

The variance of M^-'s and N^'s scale as 1/g 2 . When adding an extra electron of 
spin a in the orbital j, the variation of the residual interactions [cf. (6.7)] 

icr i 

involves in principle a summation over all occupied orbitals. The energy range within 
which fluctuations take place is however given by the Thouless energy -Erh- Thus when 
evaluating the fluctuations of 5E Rl , only a number scaling as g — E^/A of il% and 
Nij should be considered as independent. As a result the variance of 5E 111 scales as 
g x 1/g 2 = 1/g (and the r.m.s. as l/yjg). 

The contributions of the residual interactions to the second difference (6.2) then 
depend on whether or not the iV'th and N + l'th particles occupy the same orbital 
j. In the former case (assuming for simplicity that all orbitals are either doubly 
occupied or empty before the JV'th electron is added) only Mjj, survive the difference 
5 2 E m = 5E m [N^(N+l)]~5E KI [(N-l)^N]. The variance of 5 2 E K1 then scales as 1/g 2 . 
This implies in particular that the Dirac delta peak visible in figure 11 is only marginally 
affected by the fluctuations of the residual interaction. If however the last electron is 
added into a different orbital, there is no cancellation between SE m [N — * (iV+1)] and 
5E RI [(N-1) -> N}. Typical values of 5 2 E m then scales as 5E R \ i.e. as l/^/g. For 
relatively small dots (N ~ 100) this can lead to fluctuations of the order of 10% of a 
mean-level spacing. The fluctuations of the residual interactions therefore give rise to 
a contribution to the fluctuations of the second energy difference 5 2 E N which, for iV 
even (and more precisely when the iV'th and N + l'th electrons are added to different 
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orbitals), is parametrically larger (r.m.s. ~ l/y/g) than the fluctuations due to higher- 
order terms, and numerically larger than the scrambling and gate effects discussed in 
section 6.3.1. 

6.3.3. Finite temperature Up to this point, we have considered scrambling and gate 
effects, which can be considered as due to variations of the one-particle part of the 
universal Hamiltonian, and fluctuations of the residual interactions terms which are 
genuine corrections to this Hamiltonian. They give rise to fluctuations with a r.m.s. 
scaling§ as 1/ y/g and prefactors which, for experimentally realistic parameters, are 
< 1%A for scrambling, about 7% A for gate effects, and of the order of 10% A for 
the residual-interaction fluctuations. 

These figures have to be compared with the r.m.s. of the experimental noise which, 
for the dots studied in Patel et al. (1998) is equal to 8%A for dot 1 (whose data have 
been used to construct figure 7 and 8) 17%A for dot 5, 22%A for dot 4, and 40%A or 
above for the four remaining ones. 

For dots with a reasonably low level of experimental noise, these extra fluctuating 
terms (namely scrambling, gate effects, residual interaction fluctuations, and noise itself) 
can blur the sharp features of distributions such as the one depicted in figure 11. 
Nevertheless, they cannot prevent the distribution of the second ground-state energy 
difference S 2 E^ to remain bimodal, and therefore incompatible with the experimental 
peak-spacing distributions. 

Usaj and Baranger pointed out that an analysis done at zero temperature, and 
thus concerning only ground-state properties, was however inadequate to interpret the 
present experiments. For instance, the dots studied in Patel et al. (1998) were at a 
temperature ranging from 20% to 60% of a mean-level spacing. If one had in mind 
that the mean-level spacing gives the scale of the first many-body excitation, such a 
difference between T and A seemed enough to justify a zero temperature approach, at 
least for the smallest of the dots (largest A). This is especially true for chaotic system, 
since level repulsion ensure that the proportion of small spacings is small (~ 8. 10 -3 
for spacings smaller than 0.2 A in the Gaussian unitary ensemble). Therefore, within 
the constant-interaction model, the probability for an excited state to be populated 
because of thermal fluctuations could be considered as negligible for the smallest dots, 
and reasonably small for the larger ones. 

The presence of the exchange term (6.22) makes however a drastic difference in this 
respect. Consider for instance 5E s= i tS=0 , the difference of energy between the S = 1 and 
the S = lowest-energy states (assuming N even). It equals on average A — 2J S , that 
is, for r s = 1 (J s ~ 0.34), about thirty percent of the mean-level spacing. Moreover level 
repulsion does not help anymore since it affects only the small spacings. If one assumes 
GUE fluctuations for the one-particle energies, the probability that \5Es=i,s=o\ < 0.2A 
is ~ 0.29. In other words, even for the smallest of the dot studied in Patel et al. (1998), 



§ Note that for scrambling and gate effects, this scaling has been properly derived only in the case of 
billiards systems, and not for the more realistic case of smooth confinement. 
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there is almost one chance out of three that both the S = and the S = 1 states are 
both significantly occupied. 

In addition, even when most of the conduction is provided by the ground-state, the 
fact that different spin states have different degeneracies, and thus different entropies 
S = \n(2S +1), will modify the conduction peak positions when T ^ 0. 

The finite temperature linear conductance near a N — 1 — > N transition 
can be obtained in the rate equation approximation as (Beenakker 1991, Meir & 
Wingreen 1992, Usaj & Baranger 2001, Usaj & Baranger 2002) 

G(y.)-Sr p ^Ms w " (639) 

OL 

where is the equilibrium probability that the quantum dot contains N electrons, 

T ( Ft) 

T a is the partial width of the single-particle level a due to tunneling to the left 
(right) lead, and w a is a weight factor given by 

»• = k^ici**- 1 )! 2 rn^m^m ■ (6 ' 40) 

In (6.40), H,QY)\tyj f ) = Ej\tyj f ), so that "j" labels the many-body states of the quantum 
dot, and F eq (j|iV) is the conditional probability that the eigenstate j is occupied given 
that the quantum dot contains N electrons. 

From this, one can derive the temperature dependence of the peak positions 
(maxima of G(V g )) when the sole ground states have a notable occupation probability 
as (Usaj & Baranger 2001, Usaj & Baranger 2002) 

2Sg + 1 



/i = Eg - E*~ l -k B T In 



+ 2(C g E c /e) 5V g (6.41) 



which, up to a factor 1/2, amounts to replacing the ground-state energies Eq by the 
free energies F N = Eg + k B TS. 

The general case (conduction through more than one state) leads to more 
complicated expressions, but can be obtained (Usaj & Baranger 2001, Usaj & 
Baranger 2002). Both effects - conduction through excited states and displacement 
of the peaks due to the entropy - affect the spacing distribution: at k B T ~ 0.2A they 
already dominate the fluctuation of the residual interactions. This is illustrated in 
figure 12 - taken from Usaj & Baranger (2001) - where it is seen that, once scrambling, 
gate effects, fluctuation of the residual interactions, are already taken into account, 
raising the temperature from 10% to 20% of the mean-level spacing suppresses the bi- 
modality of the peak-spacing distribution. Already at T = 20%A, this distribution is 
in reasonable agreement with the experimental result shown in figure 7. As soon as the 
temperature is further increased (e.g. T = 30%A for the data in figure 8) or the level 
of experimental noise becomes significant, one should expect little (or none) even/odd 
effect or asymmetry in the peak-spacing distributions. 
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Figure 12. Finite temperature Coulomb-blockade peak-spacing distribution obtained 
once, in addition to the universal Hamiltonian, scrambling, gate effects, fluctuations of 
the residual interactions, and, more importantly, finite temperature effects, are taken 
into account. The dotted line correspond to N odd, the dashed line to N even, and 
the solid one to the total (odd plus even) distribution. N = 500 (g w 6), Js — 0.32A 
and ksT = 0.1 A (0.2A) in the top (bottom) plot. This figure is taken from Usaj & 
Baranger (2001). 



6.3.4- Non-chaotic dots As we just have seen, starting from a universal-Hamiltonian 
description, the peak-spacing distributions measured by Patel et al. can be understood 
within an approach where finite temperature effects and (experimental) noise - except 
for the quietest dot - are the most important ingredients. Conversely, the fact that noise 
and temperature play an important role implies that other features of the modelling are 
less critical, and in that sense are not really probed by these experiments. 

The size of the dot studied by Luscher et al. (2001) is however significantly smaller, 
leading to a ratio T/A ~ 5% such that temperature cannot be responsible for the 
absence of bi-modality in the peak-spacing distribution. The noise level in Luscher 
et al. (2001) is not quoted, but is presumably small enough not to affect significantly 
the spacing distributions either (Ensslin 2007). The approach based on the universal 
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Hamiltonian plus corrections does not therefore appear to be able to explain the 
distributions shown in figure 9. Other mechanisms have to be introduced to understand 
the statistical properties of peak-spacings. 

Interestingly, one issue, to which usually little attention is paid, is to know whether 
using a chaotic model for the dynamics of the electrons in the dot is adequate. 

There are in fact good reasons to focus on chaotic dynamics, the main one being 
that the behaviour of chaotic systems is universal (in the sense that this behaviour does 
not depend on the details of the system as long as this latter is in the chaotic regime). 
Working in the chaotic regime, both experimentally and theoretically, one (happily) 
avoids messy system-specific considerations. 

It remains nevertheless that, even when some care has been exercised, as in Patel 
et al. (1998), to bring the quantum dots in the chaotic regime, it is extremely difficult to 
design a system with smooth confining potential that is without question in the chaotic 
regime: No known smooth Hamiltonian system is mathematically proved to be in the 
chaotic regime; Furthermore, as any numerical simulations of the classical dynamics of 
low dimensional systems readily demonstrates, it is extremely difficult, especially in the 
presence of magnetic field, to choose a confining potential without any obvious regions 
of regular motion. 

Clearly the point is not to decide whether any experimental system in the chaotic 
regime in the strict mathematical sense, but, for a given setup, how close the system is 
from fully developed chaos, and, from a general perspective, how much the predictions 
based on the assumption that the dynamics is chaotic are robust with respect to that 
hypothesis. In other words, the question which has to be answered is how much chaotic 
dynamics are representative, for what concerns the residual interactions, of the larger 
class of mixed (i.e. partly chaotic partly regular) dynamics that one is going to find in 
practice in experiments. 

One consideration may indicate that residual interactions, and as a consequence 
the ground-state spin and peak-spacing distributions, are sensitive to the degree of 
regularity of the system under investigation. Indeed, using for simplicity the short 
range approximation V^ortrange^ — r ') — ^o^l 7 " — r ')i one can express the completely 
diagonal matrix element as Van = Ma = v^ 1 f dr\Lpi{r)\ A . Up to the factor u^ 1 , Ma 
is therefore the inverse participation ratio of the state ipi, a measure of how much this 
latter is localised. 

Of all possible dynamics, chaotic systems are however the one showing the 
least localisation, as their wave-functions are spread out uniformly across the whole 
energetically accessible domain. Integrable and mixed systems on the other hand 
display various forms of phase space localisation (in the sense that the Wigner or Husimi 
transforms of the wave-functions are concentrated in some portion of the phase space). 
The mechanisms underlying this kind of localisation (not to be confused with Anderson 
localisation) range from the relatively trivial quantization on invariant torus to more 
subtle effects of partial barrier ||, but are known in any case to be quite pervasive in 

|| For instance cantori (MacKay et al. 19846, MacKay et al. 1984a, MacKay et al. 1987) or stable and 
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mixed systems, even when the phase-space proportion of genuinely regular motion is 
not very large. 

What this considerations imply is that although chaotic systems represent most 
presumably a universality class as far as residual interaction effects are concerned, they 
might be the class for which residual interactions are the least effective when considered 
from the more general point of view of the range of possible dynamics. To go beyond 
the kind of qualitative arguments used above is made complicated by the lack of simple 
models (such as the random-plane-wave model) for the wave function fluctuations of 
mixed systems. What is possible however is to consider some specific example of smooth 
Hamiltonian systems showing, as a function of some external parameter, various kinds of 
dynamics. One then can see in these particular cases whether the actual wave-functions 
fluctuations used as an entry for the computation of the residual interactions induce 
strong modifications with respect to the chaotic predictions for the ground-state spin or 
addition-energy distributions. 

A convenient choice, introduced by Ullmo et al. (2003), is the coupled quartic 
oscillator system (r = (x,y), r = \r\)) 



x 



4 



H = ± + a(A) (^y + by 4 + 2Xx 2 y 2 j . (6.42) 

Here, b = 7r/4, a(A) is a convenient scaling factor chosen so that the mean number 
of states with energy smaller than e is given by N(e) = e 3 / 2 and A is the coupling 
between the oscillators. The parameter k is chosen such that time-reversal invariance 
is completely broken. 



Table 1. Probabilities P(S = 2) , P(S = 5/2) to find a spin two (even N) or five halves 
(odd N) ground-state, and average value (SS) of the ground-state spin augmentation 
(5S = S or (5— 1/2) for even or odd number of particles, respectively), for the various 
dynamical regimes (values of A) with k = 1.0 and J s — 0.4. The last column is the 
RMT/RPW prediction (Ullmo & Baranger 2001) valid for the hard-chaos regime. This 
table is taken from Ullmo et al. (2003). 



A 


+0.20 


-0.20 


-0.80 


RMT/RPW 


P(S = 2) 


0.13 


0.16 


0.07 


0.01 


P(S = 5/2) 


0.08 


0.10 


0.02 


0.00 


(SS) 


0.51 


0.54 


0.38 


0.23 



In figure 13, values of sets of diagonal terms are represented for one symmetry 
class in various dynamical regimes, showing as expected that their behaviour is very 
sensitive to the nature of the classical dynamics. Table 1 further shows that this 
also affects drastically the spin distribution. Indeed, for the model (6.42) there is, 
in the mixed regime, a significant proportion of ground state spin 2 or 5/2, whereas, 

unstable manifolds (Bohigas et al. 1990, Bohigas et al. 1993)). 
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Figure 13. Inverse participation ratio as a function of the orbital index for the 
symmetry class (+,+). From top down: (A, n) — (+0.20,1.00) [nearly integrable], 
(-0.20,1.00) [mixed], and (-0.80,1.00) [mostly chaotic]. This figure is taken from 
Ullmo et al. (2003). 



as predicted by the universal-Hamiltonian approach, such "large" spins are essentially 
absent in the chaotic case. As discussed in Ullmo et al. (2003), the presence of a 
large proportion of non-trivial spin is associated with distributions of spacings 5 2 En 
which differ significantly from the ones expected in the hard-chaos regime. This 
example illustrates that the question of non-chaoticity, although it has received less 
attention than other issues such as scrambling or fluctuations of residual interactions, 
may play a significant role in explaining the difference between universal-Hamiltonian 
predictions and experimental observations, and possibly be as important as the question 
of temperature. 

7. Conclusion 

In this review, I have introduced some of the tools from the field of quantum chaos which 
may be applied to the understanding of various many-body effects in mesoscopic physics. 
Among these tools, a first group is related to the semiclassical approximations of the 
Green's function. They provide a link between the quantum properties of fully coherent 
systems and the classical propagation of trajectories within their classical counterparts. 
This link gives a very intuitive picture for many quantal properties of interest. In 
particular it makes it possible to introduce naturally the classical probability P^(r, r', t) 
(cf. the M-formula (2.24)) which is central in the understanding of diffusive or chaotic 
quantum dots. Another set of tools is related to random-matrix theory, and the closely 
related random-plane-wave models, allowing for a statistical description of individual 
eigenstates of classically chaotic systems. 
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The main body of the review, however, has been devoted to the application of theses 
tools to a selected set of examples of physical interest. In particular, I have discussed 
in detail how the "semiclassical-Green's-function-based" approximations can be used to 
compute the interaction contribution to the orbital magnetic response of mesoscopic 
systems. For diffusive dynamics, I have shown that it was possible to recover in this 
way, in a simple and transparent manner, most of the results derived previously by 
diagrammatic perturbation techniques. This approach made it possible furthermore to 
address the ballistic regime, and in particular non-chaotic geometries which are beyond 
the scope of applicability of the more traditional methods. Another illustration was 
provided by the mesoscopic fluctuations of Kondo properties for a magnetic impurity in a 
bounded, fully coherent electron gas, in the temperature regime such that a perturbative 
renormalization group treatment can be applied. Here again, semiclassical techniques 
made it possible to relate some relatively non-trivial properties of the quantum system 
to the propagation of classical trajectories. 

In the last section, turning to an energy scale ~ A for which the semiclassical 
approximations of the Green's function are generally not applicable, I have discussed how 
the statistical properties of individual wave-functions derived from random ensembles 
could be applied to the study of peak-spacings and ground-state-spin statistics in 
Coulomb-blockade experiments. In the chaotic regime, I have presented in particular 
the "universal-Hamiltonian" description derived by Aleiner and coworkers, in a leading 
semiclassical (i.e. here 1/g) approximation, from a generic random-matrix argument. I 
have furthermore shown how the parameters of this universal Hamiltonian, as well as 
the evaluation of the leading corrections, can be obtained from the random-plane-wave 
models. Finally I have discussed some important modifications expected for non-chaotic 
dots, and their possible experimental relevance. 

From these various examples emerges a general picture of what the "quantum chaos" 
based approach advocated in this review brings to the understanding of the many-body 
physics of mesoscopic systems. I would like, for the remaining of this concluding section, 
to discuss it on more general grounds. 

In this respect, one important class of systems is such that either the dynamics is 
genuinely diffusive (i.e. disordered systems), or it is assumed that there is no interesting 
difference between the actual dynamics and (some limit of) a diffusive one. Ballistic 
chaotic systems under some circumstances fall into the latter category. Furthermore, 
extension of the nonlinear super-symmetric cr-model to billiard systems with diffusive 
boundaries permits addressing many important statistical properties of chaotic systems 
(Blanter et al. 1998, Blanter et al. 2001 b) (as long as the former represents a good 
model for the latter). In these cases, compared to traditional diagrammatic (Altshuler & 
Aronov 1985, Akkermans & Montambaux 2007) or more modern (Efetov 1999, Blanter 
et al. 1998, Blanter et al. 2001 b) techniques, the approach proposed in this review 
provides mainly an alternative route to the results. Its main advantage is then that 
this alternative route is closer to the actual physics considered. Indeed, most of the 
relevant parameters (Thouless energy, transport mean free path), and more generally 
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Figure 14. Illustration of an "Hikami box" (shaded region in this figure). In this 
particular example the two trajectories switch from following the same path to following 
time reversal path. Such configuration would typically be found in the calculation of 
weak localization. 

the physics (in the sense of extracting scales and qualitative behaviours) of these systems 
are usually discussed in terms of classical quantities (classical probability of return, time 
to diffuse to the boundary, etc.). The semiclassical approach developed here makes it 
therefore possible to perform the calculations with the same "language" as the physical 
discussion, and is therefore usually more transparent and intuitive. 

For disordered systems this advantage has however to be balanced with the 
respective strengths of the other approaches. The super-symmetric nonlinear cr-model 
techniques (Efetov 1999) for instance can address stronger disorders, for which Anderson 
localisation begins to develop, whereas semiclassical methods, as well as as diagrammatic 
perturbations, are limited to the "good conductor," i.e. diffusive, regime. In the same 
way, when the properties of individual wave functions are considered, and more generally 
when the physics of interest takes place on the scale of the mean level spacing A, 
super-symmetric nonlinear cr-models provide more systematic (and controlled) ways of 
performing calculations in the diffusive regime. The manipulation of the random-plane- 
wave models may on the other hand sometime involve a little bit of artwork, but is 
more flexible in nature and can be adapted to other dynamical regimes (Backer & 
Schubert 2002a, Backer & Schubert 2002 b). 

Another issue which has to be kept in mind when addressing disordered systems is 
the one of Hikami boxes (Hikami 1981). In a semiclassical language (see section C3.1 of 
Akkermans & Montambaux (2007) for a more detailed discussion) , a Hikami Box can be 
described as a restricted region of space where two pairs of trajectories following closely 
each other change partner. An illustration is given in figure 14. 

Because the length scale associated to Hikami boxes (namely the mean free path) is 
well separated from the classical scales associated to diffusion, their computation are not 
particularly difficult in the diagrammatic approach. Hikami boxes turn out to be more 
delicate to handle in a purely semiclassical description, in particular because questions 
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of current conservation have to be treated with care. Recent works $ show that these 
difficulties can actually be overcome, but the level of technicality involved becomes then 
at least as high as with the diagrammatic approach. The usefulness of the semiclassical 
approaches when processes related to Hikami boxes are involved is then not anymore 
their simplicity but the mere fact that they can be applied to a wider class of problems, 
being not by construction limited to disordered problems. 

For questions related to many-body effects in mesoscopic systems, the quantum 
chaos based approach described in this review serves therefore two purposes. One is 
to provide a way to perform needed calculations which remains very closely related to 
the qualitative physics discussed, and is therefore very transparent. In this respect it is 
presumably the best first step into this field for either theoreticians or experimentalists. 

It is, furthermore, the only way to address any experimental system which cannot 
be properly described by diffusive or fully chaotic dynamics. From this perspective, it 
is interesting to come back to the interpretation of the experimental results of Patel 
et al. (1998) and of Luscher et al. (2001). In the first case (Patel et al. 1998), some 
care has been exercised to bring the dots under study in the chaotic regime (although 
it is presumably not genuinely there). Furthermore, the value of r s (i.e. the strength 
of the interactions), and the relatively high (> 20% A) temperature make it such that 
experimental data are compatible with a chaotic modelling of the dynamics (Usaj & 
Baranger 2001, Usaj & Baranger 2002). On the other hand, the dots studied by 
Luscher et al. have an essentially square shape, and it is much less likely that the 
associated dynamics are in the fully chaotic regime. Furthermore a smaller value of r s , 
and more importantly of the temperature (~ 0.05A) make their data incompatible with 
a chaotic modelling even if one includes the effects of the gate, of scrambling, of residual 
interaction fluctuations, and of temperature. 

As seen in the last subsection of Sec. 6, the qualitative behaviour of non- 
chaotic systems can, in some circumstances, be drastically different. The experimental 
results of Luscher et al. (2001) provides at least one unambiguous example where 
assuming diffusive or chaotic behaviour is not a reasonable starting point. It therefore 
demonstrates the interest of having an approach not making too much of an assumption 
concerning the nature of the classical dynamics within the system, so that at least 
the question of whether anything new, or interesting, could happen in other dynamical 
regimes could, at least, be asked. My hope is that this review provides a step in this 
direction. 

\ Sec for instance (Alcincr & Larkin 1996, Richtcr & Sieber 2002, Adagideli 2003, Muller et al. 
2004, Stefan Hcuslcr & Haake 2004, Muller et al. 2005, Whitney & Jacquod 2005, Rahav & 
Brouwcr 2005, Heusler et al. 2006, Jacquod & Whitney 2006, Whitney & Jacquod 2006, Brouwer 
& Rahav 2006, Muller et al. 2007). 
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Appendix A. Screening of the coulomb interaction in a generic quantum 



The purpose of this appendix is to address the question of screening of the Coulomb 
interaction in a generic quantum dot. This issue has been studied in the case of diffusive 
systems (Blanter et al. 1997, Aleiner et al. 2002), and I will essentially follow here the 
approach proposed in section 2.3.2 of Aleiner et al. (2002). The main motivation will 
be to illustrate how the semiclassical approach presented in section 2 can be used to 
generalise some results derived in the more traditional framework of diffusive systems. 
Here however the exercise will turn out to be somewhat academic since I will need to 
perform at one point an "uncontrolled" (read incorrect here) approximation, leading in 
the end to an unphysical result (the equations (A.31)-(A.32) in place of the expected 
ones (2.51)-(2.52)). Nevertheless, the calculation is in itself instructive enough to be 
detailed. Identifying clearly the uncontrolled step of the derivation may furthermore 
helping clarifying the condition of applicability of the original result derived for the 
diffusive regime, which might turn out to be useful. 

In the bulk, the RPA screened interaction, is obtained by considering the Dyson 
equation for the dressed interaction (see the discussion in section 9 of Fetter & Walecka 



V AxoBSoA {r u r 2 ,u) = V conl {r 1 -r 2 )+ / dr / dr'V ^(r 1 -r)U(r,r\u)V dreBf ^(r',r 2 ,u;) ,(A.l) 



which is exact if all the one-particle irreducible diagrams are included for the polarisation 
operator n but gives the RPA approximation if only the (lowest order) bubble diagram 



is kept. G(r,r',uj) = <d(uj)G R (r,r',uj) + 0(— cu)G A (r, r', oo) is the unperturbed time 
ordered Green's function, with 6(x) the Heaviside function, and g s = 2 is the spin 
degeneracy factor. In the zero frequency low momentum limit one gets (in the bulk) 
n°(r,r', uo — 0) ~ — g s u 5(r — r'), with u the local density of states per spin (2.13). 
Inserting this expression for n° in (A.l) gives (2.46)-(2.47). 

Let us consider now a mesoscopic systems, and assume that its typical dimensions 
are much larger than the screening length. One then expects that the residual screened 



dot 



(1971)) 





(A.2) 
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Coulomb interaction should be very similar to the one in the bulk, and it is therefore 
natural to approach the question from the same viewpoint. In that case however 
the Green's function are not known exactly, so one needs to resort to semiclassical 
approximations of G R,A in the expression of II . The difficulty encountered then is that 
semiclassical approximations are valid for high energies (high a;), and in particular one 
cannot expect the semiclassical expressions for G(r, r', uo) to be accurate if uj is not 
much larger than the mean level spacing A of the system. 

Following Aleiner et al. (2002), the idea is then, in the spirit of the renormalization 
group approach, to integrate out only the "fast variable" (high-energy part) for which 
a semiclassical approximations can be used, and to deal with the low energy physics 
by some other methods (based for instance on a random-matrix description (Murthy & 
Mathur 2002, Murthy & Shankar 2003)). Using the exact expression for the polarisation 
bubble 

If (r, r>, u,)=g,J2 e(- C „ £ „.) y " (r> f toM^ (- S gn( £ „)) (A.3) 

nn' 

with (e n ,<f n (r)) the one-particle energies and eigenstates, we see that this can be 
achieved by restricting the sum in the above expression to pair (n, n') such that at least 
one energy is outside a band centered at the Fermi energy ep (taken as the origin of 
energies) and of width e* chosen such that A <C e* <C E T h, and which precise value (once 
in this range) is expected to be irrelevant. Up to an unimportant boundary term, this 
is equivalent to restricting the sum to particle-hole energies e n / — e n larger (in absolute 
value) than e*. Introducing Yi R,A {r, r', uj) == lim r) _ >0 + n°(r, r', uj ± irj) the retarded and 
advance polarisation bubbles, one can therefore write the polarisation operator in which 
only the fast modes are integrated out as 

n e .(r / ,r,fi> = 0) = ^ F / ^[U R (r,r',u;)-n A (r,r',u;)]e(\u;\-e*).(AA) 

The insertion of fl e * in (A.l) will then give the effective interaction describing the low 
energy (< e*) physics of the quantum dot. 

Appendix A.l. Calculation of the polarisation loop 

Let us first consider positive energies uj > 0. Noting that phase cancellation is possible 
only for the product G A G R , but not for G R G R or G A G A , one has 

/o j w / 

— G R (r, r>, uj' + uj + irj)G A (r', r, uj 1 ) (A.5) 

and 

U A (r,r',uj) = U R (r',r,ujy . (A.6) 

Using the semiclassical expressions (2.3) and keeping only the diagonal approximation 
in which a trajectory j is paired with itself to cancel the oscillating phases, one gets 

Air 2 

[G R {r,r\uj+J)G A {r\r,J)] A ^ = £ T | Dtf exp + uj') - S 3 {J))/h] .(A.7) 
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In this equation, one would like then to perform a Taylor expansion of the action 

(Sj(u + J) - Sj(u/)) = {dSj/de)uj = tju (A.8) 
where the last equality comes from (2.7). Inserting the unity J °° 5(t — tj) we obtain 

a _2 roo 

[G R (r, r', uj + uj')G a (t', r, uj')] diag = ^ J dt £ \Dj\ 2 5(t - tj) exp [itu MA.9) 

2vrz/ (r / ) 



/ dtP^(r,r',t)exp\itcu/h] (A.IO) 
Jo 



h 

= 2nu (r)P^r,r',uj) , (A.ll) 

where we have used the M-formula (2.24) and is the Fourier transform of the classical 
probability P^. Interestingly enough [G r (uj + uj')G A (uj')]dia,g is independent of uj', so that 
finally 

U R (r,r',uj) = -iug s v {r')Pd{r,r',u) . (A.12) 

Note that the fact that we have computed 11^, i.e. that u = uj + irj, is what is 
making the Fourier transform in (A.IO) convergent. If we had computed U A the above 
approach would have lead to divergences. I1 A should therefore be derived from 11^ using 
(A. 6), giving 

Tl A (r,r',uj) = tcogMr)P^r',r,uj) . (A.13) 

For negative u, Ii A should be calculated first and 11^ derived from it with (A. 6), leading 
to the same result. 

Here, one rather important remark is in order. The expression (A.8) assumes 
obviously that uj is small. This is usually not a significant constraint since the actions 
Si are classical quantities, so that the relevant scale is the Fermi energy (or bandwidth) 
6f- It is therefore enough that uj ep to apply (A.8). However the integral in the left 
hand side of (A. 4) is not limited to the neighbourhood of the Fermi surface. Replacing 
Y[R,a ky f. ne approximations (A.12)-(A.13) will be incorrect on the edge of the energy 
band, which will be associated to short distances \r — r'\ < A F . This will be the cause 
of the problems we shall encounter later on. Let us ignore this issue for the time being, 
and come back to this discussion when it will become obvious that the results obtained 
in this way are unphysical. 

Then, inserting (A.12)-(A.13) into (A. 4) and writing O(x) = 1 — Q(—x) we get 

— [u (r')P^r,r',uj) + u (r)P: i (r,r',uj)\ (1 - 6(e*- | w |)) .(A.14) 

The term proportional to one in the integrand of (A.14) gives rise to 
j(duj/27r)P^(r,r',uj) = P^(r,r',t^= 0). To evaluate the remaining term, it is use- 
ful to discuss the weight function @(e* — |a>|). Its precise form is irrelevant here, and, 
rather than the actual Heaviside step function, I shall assume that @(e*— \uj\) is actually 
a smooth function Q e *(uj) which is zero for \uj\ 3> e* and one for <^ e*. To fix the idea 
one can think for instance of Q e *(uj) = exp(— (l/2)(c<j/e*) 2 ), but this precise form will 
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not play any particular role. If one introduces Q € *(t) the Fourier transform of Q t *(u), 
one has, with t* = H/e* 

Q e *(t) ~ 1/t* forf«t* (A. 15) 

= for t>t* (A. 16) 

dt e e * (t) = e e * (co = o) = i . (a. 17) 



o 



Assuming furthermore that e *(i) is a positive function (this hypothesis can be easily 
relaxed), we see that B e *(t) is a density probability (since it is positive and normalised 
to one) which selects trajectory shorter than t*. We thus can write 



n e .(r',r,a) = 0) = -g s 



v Q {r')P^ry,t = Q)-\{Mr'){Pl x {r,r^ ,(A.18) 



where (f{t))t<t* == / °° dt f(t)Q e *(t) is the average over time t lesser than t* of the 
function f(t). 

Now P^(r,r',t = 0) = S(r — r'). Furthermore, the condition A < e* < E Th 
is equivalent to < f < tn, with tjj = h/A the Heisenberg time and tf the 
time of flight across the system (for ballistic systems) or time needed to diffuse to 
the boundary (for diffusive systems). We see that the choice of e* is made precisely so 
that i) semiclassical approximations are valid, but also ii) that most of the range [0, t*} is 
such that for diffusive or chaotic systems (the case of integrable or mixed system should 
be investigated in this respect), the motion can be assumed randomised. Assuming 
ergodicity we can therefore write 

This eventually leads to 

n e *(r',r,^ = 0) = -g s [u (r)5(r - r') - u (r)u (r')A] , (A.20) 

where one recognise the first term as the zero frequency low momentum bulk 
polarisationn° ulk (r', r, Cu = 0) = — g s z/ (r)<5(r — r'), and I will denote by 

n i ," W (^o(r')A (A.21) 

the remaining long range part. For billiard systems for which vq{t) = („4A) _1 = const., 
with A the volume of the system, (A.20) is for instance exactly the equation (60) of 
Aleiner et al. (2002). 

Appendix A. 2. Self- consistent equation 

In the bulk, both the Coulomb interaction (2.42) and the polarisation operator n° ulk 
are translation invariant and the Dyson equation (A.l) can be solved in momentum 
representation as 

1>M<<7) = i rrl ? ; • < A ' 22 » 

i - Kouitgjng^g) 
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The resulting interaction is then short range, effectively much weaker than the original 
Coulomb interaction, and is therefore well adapted for a perturbative treatment. 

The difficulty one encounters in the mesoscopic case is twofold. First, lack of 
translational invariance for makes in principle (A.l) impossible to be solved in 
closed form for a generic spatial variations of vq{t). Second, we know that even at the 
level of electrostatics, the effects of the interactions cannot be small since they will at 
minima rearrange considerably the static charges within the system. Therefore, even 
if (A.l) could be solved, there is little chance that the resulting dressed interaction 
could be effectively used in a perturbative approach starting from the non-interacting 
electrons Hamiltonian. 

For a quantum dot with a fixed number (N + f ) of electrons, one way to solve both 
of these difficulties is to derive a self-consistent equation following one of the standard 
derivation of the Hartree Fock approximation (Thouless 1961). For this purpose, let 
us note that any one-body potential U(r), can be written formally as the two-body 
interaction 

V(r,r') = ^(U(r) + U(r')) (A.23) 
since, using for instance a second quantization formalism 

- J drdr'^(r)i/\r')V(r,r')if(r')if(r) = j ' dr&(r)U(r)V(r) . (A.24) 

As a consequence, the total Hamiltonian, as well as the formalism presented in the first 
part of this appendix, are unmodified if the confining potential U ext (r) and the Coulomb 
potential V cou \(r,r') are respectively replaced by 

U(r) =U cxt (r) + U { r) (A.25) 
V(r, r') = Koui(r, r') - V(r, r') (A.26) 

One can then now use the freedom in the choice of the function U (r) to simplify the 
Dyson equation. In particular, if we can impose that 

J drj dr'V(r 1: r)n Lr .(r, r')V dlcsscd (r' , r 2 ) = , (A.27) 

the Dyson equation (A.l) would then have the usual "bulk-like" form 

Kiressed(n, r 2 ) = V(r u r 2 ) - J dr J dr'V( ri , r)n° ulk (r - r')V diessed {r' , r 2 ) , (A.28) 

which, if Vo{r) and U(r) vary slowly on the scale of the bulk screening length k~ 1 has 
the same solution (2.46)-(2.47) as in the bulk. 

Now, equation (A.27) might seem at first sight difficult to solve since it involves the 
unknown function Vdresscd(^ / , 7*2)- However, since U\_ T ,(r,r') does actually not correlate 
r and r', the two integrals in (A.27) actually decouple, and a sufficient condition to 
solve this equation is that J drV(r\ — r)v{r) = 0, i.e. 

J drV coul ( ri - r)u(r) = ^J drv{r){U{r) + U(r x )) . (A.29) 
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The constant term J dru(r)U(r)/N is irrelevant here as it can be eliminated by a 
constant shift of U . One therefore obtain in this way the self-consistent equation 

U(n) = NA J drV coul (n - r)u(r) . (A.30) 

In other words what we have obtained for the self-consistent potential are the equations 

U mf (r) = U cxt (r) + NA J dr'v (r')V coul (r, O (A.31) 

Mr) = J J^-A^ - U ^ r ) ~ P 2 / 2m ) ■ ( A - 32 ) 

Appendix A. 3. Discussion (C&P) 

The above result has a few nice characteristics, and one, rather unpleasant, feature. On 
the bright side, we see that it allows to clearly separate the two consequences of the 
long range Coulomb interaction: the appearance of a self-consistent potential (obtained 
from equation (A.30)) on the one hand, and on the other hand the screening process 
leading to the usual "bulk" form ((2.46)-(2.47)) of the screened interaction (provided 
that U(r) (and thus V(r, r')) varies slowly on the scale of the screening length). 

What makes (A.30) less useful however is that it is obviously incorrect. Indeed we 
know that whatever self-consistent equation we end up writing, it should contain in some 
approximation the electrostatic equilibrium of the problem. This is not the case here. 
If the self-consistent potential U m {(r) = f U ext (r) + U(r) obtained from (A.30) is well 
approximated by a constant (giving for instance a billiard system with weak disorder 
as was considered in Blanter et al. (1997) and Aleiner et al. (2002)), and assuming 
(N + 1) ^> 1, one can do the replacement NAu(r) — > n(r) in (A.30) and write instead 

U( ri ) = J drV coul ( ri - r)n(r) , (A.33) 

i.e. (2.51)-(2.52), which is just the Thomas Fermi equation, from which plain 
electrostatic is obtained by neglecting the kinetic energy term TtfH in (2.48). However 
for a generic confining potential U cxt (r), solutions of (A.30) will not in general be an 
approximation of the solution of (A.33) 

What we see is that, in some sense, equation (A.30) is "aware" of the properties of 
the system near the Fermi energy (the density of states v(r)), but misses the relevant 
information at large energies, of the order of the bandwidth. This is to be expected 
since the polarisation operator II* (A. 20) involves only the local density of states at 
the Fermi energy z/(r). This can be tracked back to the approximation (A. 8) where 
the action S(u>) has been linearised near the Fermi energy, eliminating in this way any 
information relevant to the large (i.e. ~ e^) energies. 

Appendix B. Magnetisation and persistent current 



In this appendix, I re-derive briefly (for sake of completeness) the basic expressions (3.2) 
and (4.26) of the magnetisation and persistent current which are the starting points of 
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the discussion in sections 3 and 4. I follow here the presentation given in Desbois et al. 
(1998). 

Let us therefore consider a two dimensional gas (d = 2) of electrons confined by a 
potential U(r) (r = (x, y) are the coordinate inside the plan and hatz the unit vector in 
the perpendicular direction) and interacting through V(r — r'). The total Hamiltonian 
of the system is therefore expressed as H = H + H mt with 

1 



— (-«V-eA(r)) 2 + [/(r) 
2m 



#(r) 



H d £ J dr¥{r) 

the non interacting part (A(r) is the vector potential), and 

= \ J drdr'¥{r)¥{r')V{r - r')if(r')if(r) , 

the interacting part (which is however not going to play any role here). One may 
furthermore introduce the current density operator 

j(r) = e^ t (r)v^(r) 

with 

v d = —i—ihV - eA(r)) 

m 

the velocity. The variation of the Hamiltonian corresponding to a variation 5 A of the 
vector potential is then expressed as 



8H = ~ J dr \j{r)5A{r) + SA(r)j(r)} . 



Appendix B.l. Uniform perpendicular magnetic field 



Let us consider first the case where the variation 5 A corresponds to a uniform magnetic 
field . 
with 



field B = 5B z z. The equation (3.2) basically state that the magnetisation M z = f (M z ), 



M z = i J dr{r x j) ■ z , 

is the variable conjugated to SB Z . Indeed, choosing for convenience the symmetric gauge 
SB Z 

and, noting that (z x r) and v commute, we have 

r) B f 

5H = --^ J dr[j(r) ■ (z x r) + (z x r) ■ j(r)} = -5B Z M Z . 

We see that M z is indeed conjugated to 5B Z , and the variation of the grand potential 
(3.1) gives (3.2). 
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Appendix B.2. Flux line 

Let us now consider a half infinite line V = {r + au; (a > 0)} originating from r and 
directed along the unit vector u = (cos 9 , sin 9 ). One can define the current operator 
across V 

?orb/_ n N def /" Jt _ _,(»*- f ) X j(r) 



/ OI >oA) = / d|r-r | v , , JV y 



or, if no current escape to infinity so that (J orb ) has no 9 dependence 

/«>„, s I J, 

Introducing the vector potential 



r — r r 



Avoir) = — $ 



def 1 ,ix(r - r ) 



2n \r — r \ 

we see that the corresponding magnetic field is B = V x A ro = $<5(r — ro)z, and 
therefore describes a flux lines $ at ro- Following the same steps as above, we find the 
variation of the Hamiltonian associated to a variation <5$ of the flux to be 

5H = -5$ I dr¥(r) — Z * ^ ~ ^ ■ vV(r) = -5$/ orb (r ) 
J 2ix \r-r \ 2 

The orbital current J orb (ro) is thus conjugated to the flux $, which, noting / = f (J orb ), 
as before directly imply (4.26). 

Appendix C. List of symbols 

• aj M : random-plane-wave coefficient (see (2.39)). 

enclosed by the trajectory j. 

• A{r): vector potential. 

• B: magnetic field. 

• /3: 1/fcnT. 

• (3rmt- random-matrix ensemble parameter. 

• C: quantum dot total capacitance. 

• C g : gate dot capacitance. 

• cLa, c aa : creation and destruction operators. 

• x'- magnetic susceptibility. 

• Xl- Landau susceptibility. 

• Xloc : local susceptibility. 

• D: diffusion coefficient. 
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• Dj: determinant describing the stability of the trajectory j (cf. (2.5)). 

• Do: original bandwidth of the Kondo problem. 

• D e g: running bandwidth of the Kondo problem. 

• A: one-particle mean level spacing. 

• e K : one-body energy. 

• Ejy: many-body ground-state energy of a N-particle quantum dot. 

• E{ni a }: energy of the many-body state corresponding to the occupation numbers {ni a }- 

• E sm (N): smooth part of E{rii a }. 

• E R1 {rii a }: contribution of the residual interactions to E{ni a }. 

• E c : charging energy. 

• E^h'- Thouless energy. 

• £ cx t : Confinement part of the Thomas Fermi functional. 

• <?coui : Coulomb part of the Thomas Fermi functional. 

• /o : Fermi liquid parameter. 

• /(e — fj,): Fermi function. 

• / x (T7Tr-): universal function describing the susceptibility for the [bulk] Kondo problem. 

• Ttf- Thomas Fermi functional. 

• g: dimensionless conductance Eq^h/ A. 

• g s : (=2) spin degeneracy. 

• G(r,r';e): (unperturbed) Green's function. 

• G A (r,r';e): (unperturbed) advanced Green's function. 

• G R (r, r'\ e): (unperturbed) retarded Green's function. 

• G^(r, r'\ e): semiclassical contribution of the orbit j to the retarded Green's function (see 
(2.3)). 

• Q(r,r';u): (unperturbed) Matsubara Green's function. 

• Ta^: partial width of the single-particle level a. 

• H: Planck constant. 

• Jo: bare coupling constant of the Kondo Hamiltonian. 

• Jeff: renormalised coupling constant of the Kondo Hamiltonian. 

• Js- Universal-Hamiltonian coupling constant (S^ ot term). 
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• Jt- Universal-Hamiltonian coupling constant (T^T term). 

• <^rpa: RPA approximation of J5. 

• k F : Fermi wave- vector. 

• k(r): h- l v / 2m(E- U{r)). 

• k: screening wave-vector. 

• t: mean free path. 

• L: characteristic size of the system. 

• Lt- thermal length. 

• L^: coherence length. 

• Af: Fermi wavelength. 

• Ao: electron-electron bare coupling constant. 

• A (A): electron-electron running coupling constant. 

• Ao: bare cutoff. 

• A: running cutoff. 

• m e : particle (electron) mass. 

• Mij-. screened interaction matrix element (see (6.8a)). 

• (M z ): magnetisation. 

• ji: chemical potential. 

• fi mc : micro-canonical distribution (2.36). 

• N: number operator. 

• Njj-: screened interaction matrix element (see (6.86)). 

• v\ oc (r,e): local density of states (see (2.1.2)). 

• fo(r,e): Weyl (smooth) part of the local density of states (see (2.12)). 

• v osc (r,e): oscillating part of the local density of states (see (2.15)). 

• Vj{r, e): Contribution of the orbit j to the oscillating part of the local density of states. 

• Vfs(r,T): thermally averaged local density of states (see (5.10)). 

• n(r): density of particles. 

• n(r,r'): non diagonal element of the particle density matrix (see (3.5)). 

• f2: grand potential. 

• Q c : Cooper series contribution to the grand potential. 
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• P^(r,r',t): classical probability to go from r' to r in a time t at energy e. 

• P^(r,r f ,uj): Fourier transform of P^(r, r', t). 

• (puir)'- one-body eigenstate. 

• Mw ( r )P) : Wigner transform of the state (p. 

• p F : Fermi momentum. 

• P nns : nearest neighbor distribution. 

• (po- quantum of flux. 

• 4>j\ magnetic flux enclosed by the trajectory j. 

• \F^: many-body eigenstate of a quantum dot with iV-particles. 

• ty a (r), if a (r): creation and destruction operators. 

• r s : electron gas parameter. 

• R(x): xj sinh(x). 

• pit): total density of states. 

• po( e ) : Weyl (smooth) part of the total density of states. 

• S = (S x , S y , S z ): spin operator (= ^ha). 

• S to t: quantum dot total spin. 

• a = (a x ,a y ,a z ): Pauli matrices. 

• S(r,r';6<j): (free) particle-particle propagator (see (3.16)). 

• £( D )(r, r'; u): diagonal part of the particle-particle propagator. 

• Y±j(r,r';uj): contribution of the orbit j to the diagonal part of the particle-particle 
propagator. 

• tf time of travel of the orbit j. 

• ig: time of flight across a ballistic structure. 

• tT'- thermal time (see (3.8)). 

• Tk'- Kondo temperature (bulk case). 

• T^: average Kondo temperature (mesoscopic case). 

• T^[up]: realization dependent Kondo temperature (mesoscopic case). 

• Tk [vp\ {T) : realization and temperature dependent Kondo temperature (mesoscopic case). 

• 7xf : kinetic energy part of the Thomas Fermi functional. 
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• 6 (it): Heaviside function. 

• U(r): one-body potential. 

• U cx t(r): one-body external potential. 

• U m f(r): one-body mean- field potential. 

• V g : gate voltage. 

• V*: gate voltage corresponding to a conductance peak. 

• V(r — r'): two-body interaction. 

• Koui( r — r ') : bare Coulomb interaction. 

• Vsc{ r — r ') : screened Coulomb interaction. 

• Khortrange { r — r ') : short range approximation (= Uq 1 5(t — r') of the screened interaction. 

• Vijki'- screened interaction matrix element (see (6.19)). 

• Vsc(q)'- Fourier transform of the screened Coulomb interaction. 

• (Kic)f.s. : Fermi surface average of V sc (q). 

• Q: Maslov index of orbit j. 
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